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Abstract 

This paper provides a systematic description of the interplay between a specific class of 
reductions denoted as cKP rm (r, m > 1) of the primary continuum integrable system - the 
Kadomtsev-Petviashvili (KP) hierarchy and discrete multi-matrix models. The relevant inte- 
grable cKP r!m structure is a generalization of the familiar r-reduction of the full KP hierarchy 
to the SL(r) generalized KdV hierarchy cKP r .o- The important feature of cKP rjm hierarchies is 
the presence of a discrete symmetry structure generated by successive Darboux-Backlund (DB) 
transformations. This symmetry allows for expressing the relevant tau- functions as Wronskians 
within a formalism which realizes the tau- functions as DB orbits of simple initial solutions. 
In particular, it is shown that any DB orbit of a cKP r .i defines a generalized 2-dimensional 
Toda lattice structure. Furthermore, we consider the class of truncated KP hierarchies (i.e., 
those defined via Wilson-Sato dressing operator with a finite truncated pseudo-differential se- 
ries) and establish explicitly their close relationship with DB orbits of cKP r . m hierarchies. This 
construction is relevant for finding partition functions of the discrete multi-matrix models. 

The next important step involves the reformulation of the familiar non-isospectral additional 
symmetries of the full KP hierarchy so that their action on cKP rm hierarchies becomes consistent 
with the constraints of the reduction. Moreover, we show that the correct modified additional 
symmetries are compatible with the discrete DB symmetry on the cKP r m DB orbits. 

The above technical arsenal is subsequently applied to obtain complete solutions of the 
discrete multi-matrix models. The key ingredient is our identification of q-matrix models as DB 
orbits of cKP r .i integrable hierarchies where r = (p q — 1) . . . (p% — 1) with p\, . . . ,p q indicating 
the orders of the corresponding random matrix potentials. Applying the notions of additional 
symmetry structure and the technique of equivalent hierarchies turns out to be instrumental in 
implementing the string equation and finding closed expressions for the partition functions of 
the discrete multi-matrix models. As a byproduct, we obtain a representation of the r-function 
of the most general DB orbit of cKPi.i hierarchy in terms of a new generalized matrix model. 

1 Work supported in part by the U.S. Department of Energy under contract DE-FG02-84ER40173 
2 Supported in part by Bulgarian NSF grant Ph-401 



The present formalism is of direct relevance to the study of various random matrix problems 
in condensed matter physics and other related areas. In particular, we obtain a new type of 
joint distribution function with an additional attractive two-body and a genuine many-body 
potentials. 
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1 Introduction 



1.1 Motivation and Strategy 

The main purpose of this paper is to present a conceptually consistent and complete description 
of a special class of integrable hierarchies of non-linear differential equations - truncated and con- 
strained Kadomtsev-Petviashvili (KP) hierarchies, identified as specific reductions of the standard 
full-fledged KP hierarchy. Another purpose is to apply this structure towards the goal of obtaining 
closed form solutions to the discrete multi-matrix models. 

There are many equivalent formulations of the theory of integrable models and, clearly, any 
account will have to be selective and focus on one or few approaches which best fit the gen- 
eral strategy. Our presentation follows a path which exhibits the connection between integrable 
models and multi-matrix models. In our case, the strategy dictates that we start with the pseudo- 
differential operator calculus of Lax operators and then establish the links to dressing operators 
and Baker- Akhiezer (wave) functions of the underlying linear problem. This in turn, allows us to 
formulate the relevant integrable models in terms of their tau-functions. It is conceptually impor- 
tant to understand that the discrete structure of matrix models is being captured by the canonical 
Darboux-Backlund symmetry in the setting of the continuous integrable models. In this spirit, 
the tau-functions are constructed explicitly as Darboux-Backlund orbits of certain simple initial 
solutions. This part of the construction is directly relevant for finding the partition functions of 
multi-matrix models. The story is, however, not as simple as it was briefly sketched above due to 
the existence of the so called "string" equation constraining the form of corresponding partition 
functions. Accordingly, the complete understanding of the connection between multi-matrix models 
and constrained KP hierarchies requires several important additional concepts. 

First of all, we need to formulate the notion of additional symmetries for constrained KP hier- 
archies and describe its interplay with the discrete symmetry structure provided by the Darboux- 
Backlund transformations. Secondly, one discovers that, within the formalism of multi-matrix 
models, the fruitful notion of equivalent hierarchies facilitates implementation of the string equa- 
tion before the final identification with the relevant Lax structure of the constrained KP hierarchy. 
Let us now briefly describe the key ingredients of our construction. 

1.2 Integrable Systems and KP Hierarchy 

Integrable systems constitute an outstanding branch of theoretical physics and mathematics (for 
the basics, see refs.(l|, ||[ [|, ^|). The list of their applications to physics is enormous: they 
have been found essential for the description of a vast variety of fundamental non-perturbative 
phenomena ranging from the hydrodynamics of D = 2 (space-time dimensional) nonlinear soliton 
waves and planar statistical mechanics to string and membrane theories in high-energy elementary 
particle physics. Furthermore, a variety of physically interesting theories in higher space-time 
dimensions can be reformulated (under plausible assumptions keeping the essential properties of the 
dynamics) as lower-dimensional (D = 2) integrable models which in the same time possess infinite- 
dimensional symmetries (see, especially, the recent developments @ related with integrability of 
effective low-energy theory of (extended) super symmetric gauge theories). On a more formal level, 
the theory of integrable systems is a corner-stone of various basic disciplines in mathematics itself, 
such as algebraic geometry (Schottky problem), group theory (quantum groups) and topology (link 
polynomials in knot theory) . For a recent account of the latter developments we refer to ref . |8| . 

The primary integrable system, based on ^fDO (the algebra of all pseudo-differential operators 
on the circle) as symmetry algebra, is the KP integrable hierarchy [||, ^] of soliton nonlinear evolution 
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equations, which plays a role of a universal arch-type integrable system. Its universality is underlined 
by the fact that the KP hierarchy (and multi-component generalizations thereof) is intimately 
related to: the discrete generalized Toda lattice hierarchy |J (which similarly plays a role of primary 
lattice integrable system); integrable systems of point particles |l(], [|; theory of Riemann surfaces 
PI , U ; theory of representations of infinite-dimensional Lie algebras (T^| . 

Since the seminal works of Sato's school Jl3| (and also of Segal- Wilson |14|]) it is well-known 
that any solution of KP hierarchy can be expressed in terms of a single function of all evolution 
parameters - the so called T-function, satisfying fundamental bilinear Hirota equations. The space 



of all r-functions is known as the universal Sato Grassmannian (see ref . [ 15 1 for a non-technical 
introduction to Sato theory). 

In the last few years, the main interest towards KP hierarchy originates from its deep connection 



fllq , |l7| , 18] with the statistical-mechanical models of random matrices ((multi-)matrix models) 
providing non-perturbative discretized formulation of string theory (for a comprehensive review, 
see refs.Q). 

One of the most remarkable achievements in the theory of (multi-)matrix models of string theory 
is the realization that the (multi-)matrix partition functions are r-functions of certain reductions 
of the full KP hierarchy. Therefore, the problem of explicit identification of the relevant reduced 
KP hierarchy in terms of its Lax operator description turns out to be of definite interest since the 
Lax formalism provides the most succinct formulation of integrable systems. In the case of discrete 
(multi-)matrix models (i.e., without taking the double-scaling limit), the pertinent reduced KP 
hierarchies belong to the class described in that which follows. 

1.3 Constrained KP Hierarchies 

The fundamental character of the KP hierarchy shows up, among other things, in the fact that via 
various reductions schemes it provides a unifying description of a number of basic soliton equa- 
tions like (generalized) Korteweg-de-Vries (KdV), nonlinear Schrodinger (NLS, or more generally 
- the AKNS hierarchy), Yajima-Oikawa, coupled Boussinesq-type equations etc.. The latter are 
embedded in different ways in various reduced KP hierarchies comprising the class of the so called 
cKP rm constrained KP hierarchies. 

This large and important class of integrable hierarchies of non-linear differential equations, 
appearing as reduced KP hierarchies, is the principal object of the present paper. The relevant re- 
duction scheme we are going to describe here generalizes (and is less restrictive that) the well-known 
r-reduction of the full KP hierarchy to the SL(r) Gelfand-Dickey (generalized KdV) hierarchy. 

The cKP rm integrable hierarchies appear in different disguises from various parallel develop- 
ments: 



Symmetry reductions |2Cj, [2l|, |22|] of the full KP hierarchy - the "eigenfunction" [23, |2J] form 
of cKP r , m (eqs.(|23H),(p3D below). 

Abelianization, i.e., free-field realizations, in terms of finite number of fields^] of both compat- 
ible first and second KP Hamiltonian structures - "m-generalized two-boson" form [25, |2(| 
and the "ratio" (or SL{r + m,m)) form p7|, |28|, |29| of cKP r;m , respectively, generalizing 
SL(r) KdV hierarchy (eq.fl23|) and eqs.(^3|)^3^) below); (see also refs.||, 

A method of extracting continuum integrable hierarchies (by taking fields associated with one 



fixed lattice site) from the generalized Toda-like lattice hierarchies [ 31 1 underlying (multi 



3 Recall that the full KP system contains an infinite number of coefficient fields. 



4 



)matrix models - "multi-boson" KP reductions (eq.( |2.52j ) below) generalizing the familiar 
two-boson KP reduction [g2|, |3^| . 

• Purely algebraic approach in terms of the zero-curvature equations for the affine Kac- Moody 
algebras with non-standard gradations p4| (interpolating between homogeneous and principal 
gradations of the Drinfeld-Sokolov approach). 

There exist explicit Miura-like transformations connecting the above mentioned different forms of 
cKP rm hierarchies (see eqs.( 2.53|) , (|2.54|) and ( 2.56| ) below). 



1.4 Darboux-Backlund Transformations, Associated Lattice Model Structures 
and Discrete Multi-Matrix Models 

The recent surprising discoveries linking the discrete multi-matrix models to continuum integrable 
systems have their roots in the formulation of the constrained KP systems possessing extra discrete 
symmetry - Darboux-Backlund (DB) symmetry, which mimics lattice shifts between neighboring 
sites. 

Originally, the classical Darboux transformation was introduced in the context of the linear 
spectral problem for purely differential operators of Schrodinger type. Namely, using a known 
eigenfunction one can transform the initial Schrodinger operator into a new one with modified 
potential and express the spectral data of the transformed operator explicitly in terms of the 
transformation-generating eigenfunction. Since auxiliary linear spectral problems are inherently 
associated with the nonlinear evolution equations of integrable systems, Darboux transformations 
play an important role in generating multi-soliton solutions of the latter (for an extensive overview, 



see ref.[35|). 

Recently, the notion of Darboux transformations was generalized to the case of pseudo-diffe- 
rential operator formulation of integrable hierarchies - the so called Darboux-Backlund transfor- 



mations (see [21] and references therein). Especially, we will be interested in auto-Backlund trans- 



formations applied to the cKP rim hierarchies - i.e., such DB transformations which map solutions 
of a given initial cKP ri?n integrable hierarchy into other solutions of this same hierarchy. Let us 
stress that in our approach it is very important that all Lax operator eigenfunctions, generating 
the DB transformations, are not Baker- Akhiezer (wave) functions (unlike the case in ref.|2l|): they 
are rather "spectral" superpositions of the latter. 

As demonstrated below in section ||, the associated Darboux-Backlund orbits, i.e., the set 
of all solutions of the subclass of cKP r i hierarchies obtained upon successive applications of 
(auto-)Darboux-Backlund transformations on a fixed initial solution, give rise to generalized two- 
dimensional Toda-lattice-like structures which, in turn, naturally arise in the context of (multi- 
)matrix string models. This generalizes the known connection between ordinary (one-dimensional) 
Toda lattice models and Backlund transformations of refs.[36, 33, p7[]. 



The eigenfunctions and r-functions of all members of any cKP rim DB orbit can be explicitly 
expressed in Wronskian form entirely in terms of the data in the initial cKP rm hierarchy (see 
eqs.( 3.40| )-( 3.42j ) and ( |3.53D -( j3.54 ) below). There exists alternative "spectral" representation for 



the r-functions of cKP r;m DB orbits (cf. eq.( p.47|) below) which turns out to be very instrumental 
in interpreting of such r-functions as partition functions of (multi-)matrix models generalizing the 
familiar ones with standard polynomial matrix potentials. 
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1.5 Additional Symmetries of Integrable Hierarchies and the String Equation 

Any integrable hierarchy is characterized by the presence of an infinite number of mutually com- 
muting flows of the associated evolution parameters ( "times" ) - vector fields acting on the space of 
the corresponding Lax operators. They form an infinite-dimensional Abelian algebra of isospectral 
symmetries of the associated nonlinear evolution equations. However, the space of all symmetries of 
the nonlinear evolution equations turns out to be much broader |3£|. The extra symmetries, which 
are explicitly time-dependent but commute with the ordinary isospectral ones, are called "non- 
isospectral" or "additional" symmetries and, furthermore, they form a non-commutative algebra. 
A particularly efficient approach to deal with additional symmetries was developed by Orlov and 
Schulman p9[ . For systematic reviews we refer to [30, ftQ| ; see also the discussions [41] of additional 



symmetries in the context of specific integrable models (AKNS, truncated KP, and generalized 
matrix hierarchies, respectively). 

In the Orlov-Schulman approach, the additional symmetries are constructed as flows (vector 
fields) acting on the space of Lax operators of the integrable hierarchy. As proved explicitly in 
refs.p2|, fjjfl , there exists an equivalent definition of additional symmetries as vector fields acting on 
the space of r-functions (Sato Grassmannian) of the corresponding hierarchy. This latter formula- 
tion allows to provide a simple interpretation of the Virasoro (and Wi +00 ) constraints on partition 
functions of (multi-)matrix models in the double-scaling limit where the underlying integrable hier- 
archies are of the generalized SL(r) KdV type. Namely, the Virasoro constraints express invariance 
of the r-functions (i.e., the partition functions) under the Borel subalgebra of the Virasoro algebra 
of additional non-isospectral symmetries (similarly for the Wi +OG constraints). In particular, the 
so called "string" equation in (multi-)matrix models is nothing but condition for invariance of the 
partition function under the lowest additional symmetry flow. 

There is, however, a problem in applying the Orlov-Schulman formulation of additional sym- 
metries in the case of the cKP r]m hierarchies. Namely, the standard Orlov-Schulman additional 
symmetries for the full unconstrained KP hierarchy are broken by the constraints when restricting 
their action to the cKP r m Lax operators. Therefore, they need to be suitably modified in order 
for the corresponding additional symmetry flows to preserve the space of cKP rm Lax operators for 
any fixed r, m . The solution to the problem is adding to the Orlov-Schulman flows appropriate 
linear combinations of the so called "ghost" symmetry flows generated by pairs of eigenfunctions 
and adjoint eigenfunctions of the cKP rim Lax operator; see ref.|44] and section |] below. 



For Darboux-Backlund orbits of cKP r m hierarchies we have the following important property. 
The "string" equation (i.e., the condition for invariance under the lowest additional symmetry flow) 
on the initial hierarchy alone turns out to be sufficient to uniquely determine the r-functions of the 
whole cKP rm DB orbit which, thus, automatically satisfy the pertinent higher constraints (under 
higher additional symmetry flows). 

Implementing the string equation constraint in the setting of the discrete multi-matrix models 
and its corresponding linear problem taking the form of Toda-like lattice system turns out to be 
greatly simplified by application of the techniques of the equivalent hierarchy method. As shown in 



[29], the constrained two-matrix Toda lattice structure augmented by the string condition can be 
re-expressed, by a simple change of variables in the set of coupling parameters as a single set of flow 
equations for only one independent matrix. This technique is extended here to the multi-matrix 
system and the resulting simplification opens the way for identification of the constrained Toda-like 
lattice system with the corresponding continuum generalized reduced KP hierarchy with the string 
equation constraint imposed on it. 

The method of equivalent hierarchies appeared also in the setting of the Toda lattice hierarchy 
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[45] and the Generalized Kontsevich Model p6f| . 



1.6 Plan of Exposition 

In section || we start with a brief review of the basics of Sato formulation of the standard uncon- 
strained KP hierarchy. Among the covered topics are: pseudo-differential Lax operator formalism, 
Baker- Akhiezer (wave) functions, r-functions, bi-Hamiltonian structure. Then we go on with intro- 
ducing the class of reduced (constrained) cKP rjm hierarchies in various equivalent forms and provide 
explicit Miura-like transformations connecting them. Also, as a byproduct we describe series of 
explicit free-field realizations of Wi +oc and algebras - the first and the second fundamental 
KP Poisson bracket algebras. 

Section ^ is devoted to the explicit construction of Darboux-Backlund solutions for arbitrary 
cKP rm hierarchies. First, after recalling the basic definition and properties of DB transformations, 
we give the derivation of the two-dimensional Toda lattice structure corresponding to any DB orbit 
of cKP r>m hierarchies. Further, we explicitly construct, in terms of simple Wronskian determinants, 
the eigenfunctions and r-functions of these general DB orbits. In particular, we provide a "spectral" 
representation of the cKP r ^ m DB-generated T-functions which will become useful for the subsequent 
identification of the latter as "eigenvalue" integral representations of the (multi-)matrix model 
partition functions. 

In the second half of section ^ we discuss a special class of constrained KP hierarchies - the 
so called truncated KP hierarchies, defined in terms of dressing operators with finite (truncated) 
pseudo-differential series. This class includes especially the subset called generalized Burgers- 
Hopf hierarchies, constructed in terms of the DB transformations of the most simple, non-trivial 
dressing operator built of one function only. We identify the truncated KP hierarchies as particular 
members of certain cKP r ^ m DB orbits and present the explicit solutions for their BA and r-functions 
which, furthermore, are shown to contain the well-known multi-soliton KP solutions. 

In section |] we study and solve the problem of modifying the usual Orlov-Schulman additional 
symmetries in order for them to preserve the space of Lax operators of any given constrained 
cKP rm hierarchy (for fixed r,m). Explicit solution is given for the Virasoro algebra of additional 
symmetries, and the generalization of the latter construction for the full Wi +oc algebra of additional 
symmetries is indicated. 

In section |5| the techniques developed in the previous sections are applied to the treatment of the 
simplest model of random matrices - the Hermitian one-matrix model. Its underlying integr ability 
structure is identified with the above mentioned special case of truncated KP hierarchy - the 
generalized Burgers-Hopf system, subject to the additional condition of invariance under the lowest 
additional symmetry flow (the "string" equation constraint). Also, it is explicitly demonstrated 
how this last single constraint uniquely fixes the Burgers-Hopf r-function (i.e., the one-matrix 
model partition function) such that the latter automatically satisfy all higher Virasoro additional- 
symmetry constraints. Also, as a byproduct, we obtain a representation of the r-function of the 
most general DB orbit of cKPii hierarchy in terms of a new generalized matrix model (cf. eq.(|5.4|) 
below) . 

Section || is devoted to a detailed study of the integrable structure underlying the discrete (Her- 
mitian) two-matrix model. After recalling the basic properties of the corresponding generalized 
Toda-like lattice system, we employ a method similar to the method of equivalent hierarchies to map 
the latter lattice system into an equivalent, but simpler one, where part of the pertinent "string" 
equation constraints are automatically fulfilled. Furthermore, by taking quantities associated with 
a single lattice site we explicitly construct, out of the modified Toda-like lattice system, the rele- 
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vant continuum integrable hierarchy which is identified as constrained KP hierarchy of cKP rm= i 
type subject to the subsidiary condition of invariance under the lowest additional (non-isospectral) 
symmetry flow. The latter identification allows us to apply the Darboux-Backlund techniques from 
section || in order to obtain the explicit Wronskian solution for the two-matrix model partition 
function. 

In section |7] we extend the treatment from section ^ to the case of arbitrary multi-matrix models. 

The concluding section |8| contains brief discussion of some further developments, in particular, 
the relevance of the present formalism for studying random matrix problems in condensed matter 
physics. 

Some useful formulas and identities are collected in the Appendix. 

Preliminary short expositions of part of the present results were previously reported in refs.|p9|, 
23,0,H,0]. 



2 KP Hierarchy and its Reductions 
2.1 Background on KP Hierarchy 

We use the calculus of the pseudo-differential operators (see e.g. [13, to describe the KP 
hierarchy of integrable nonlinear evolution equations. In what follows the operator D is such that 
[D , /] = /' with f = df = df/dx and the generalized Leibniz rule holds: 

D n f = Ys\ n \d j f)D n - j , neZ (2.1) 



3=0 



In order to avoid confusion we shall employ the following notations: for any (pseudo-)differential 
operator A and a function /, the symbol A{f) will indicate application (action) of A on /, whereas 
the symbol Af will denote just operator product of A with the zero-order (multiplication) operator 
/• 

In this approach the main object is the pseudo-differential operator: 

r— 2 oo 

L = D r + Y,v j D j +^u iJ D- ,; (2.2) 

j=0 i=l 

called (generalized) Lax operator. The Lax equations of motion: 

^ = [lI,L] , n = l,2,... (2.3) 



describe isospectral deformations of L. In (|2.3| ) and in what follows the subscripts (±) of any 
pseudo-differential operator A = J2j denote its purely differential part (A + = J2j>o ajD^) or 
its purely pseudo-differential part (A- = J2j>i a -jD~ 3 ), respectively. Further, (t) = (t\ = x, £2, • • •) 
collectively denotes the infinite set of evolution parameters from (^^) . 

One can also represent the Lax operator ( [2 .2] ) in terms of the dressing operator: 



W = 1 + J2 w nD~ n ; L = WD r W- 1 (2.4) 



n=l 



In this framework eqs.(pT3|) are equivalent to the so called Wilson-Sato equations: 

dW — — 

d n W = — = LlW -WD n = -LLW (2.5) 
dt n 
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(in what follows we shall use frequently the short-hand notation d n = gf-). Define next the Baker- 
Akhiezer (BA) function via: 



i/>(t,\) =W(e^' X) ) = w(t,X)e^' X) ; w(t, A) = 1 + ^ w i (t)X~ i , (2.6) 

where 



i=l 



£(i,A) = 5> n A n ; i 1= x (2.7) 

n=l 

Accordingly, there is also an adjoint BA function: 

oo 

^ = W*- 1 (e-^) = w*(t,X)e-^' X) ; w*(t, A) = 1 + £ <(t)A _i (2.8) 



and one has the following linear system: 



L(^(t,A)) = A>(i,A) , cW = ^W>M)) (2.9) 
L*(^(t,A)) = A>*(t,A) , 9 n V*(t,A) = -(L*)l(^(i,A)) 



Note that eqs.( |2.3| ) for the KP hierarchy flows can be regarded as compatibility conditions for the 
system (|2~9| ). 

There exists a quite natural way of developing the KP hierarchy based on one single function 
of all evolution parameters - the so-called tau function r(t) . This approach is an alternative 
to using the Lax operator and the calculus of the pseudo-differential operators. The r-function is 
related to the BA functions (|2~^)~(^9|) via: 



T ( t ) n=0 T ( t ) 

r(t, A) = T(t + . [ f 1]) e-^> = e-*™ £ ^Ml^-X-n (2 . n) 
T \ t ) n=Q T \ t ) 

where: 

[A" 1 ] = (a- 1 , \x~\ \x~\ . . .) ; [5] = . . .) (2.12) 
and the Schur polynomials p n (t) are defined through: 

e^i^^^A^p^i^s,...) (2.13) 

71=0 

Taking into account ( |2.10| ) and ( |2.6[) , the expansion for the dressing operator (2.4) becomes: 

w = t ntkMl® D -n (2 . 14 ) 

The BA functions enter the fundamental bilinear identity 

ip(t,X)ip*(t',X)dX = (2.15) 
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which generates the entire KP hierarchy. In ( 2.15| ) f^^dX is the residue integral about oo. It is 
possible to rewrite the above identity in terms of the tau-functions obtaining: 



(2.16) 



r(t - [X- 1 ])^' + [X-^e^-^'^dX = 
Taylor-expanding ( |2.16| ) in y (t — > t — y, t' — > t + y) leads to: 

(£^(-2y)p n+1 ([ J D])e5:r^ D r • r = 



(2.17) 



with [D] = (D\, (1/2)L>2; (1/3)I?3, . . .), where Di is the Hirota derivative defined as: 

b = (d m /dsf)a(t j + SjMtj - Sj )\ s=0 



D T a 



The coefficients in front of y n in ( 2.17 ) yield: 

-D\D n -Pn+l 



([£>])) t-t = 



(2.18) 



(2.19) 



Eqs.( |2.17| ) and (2.19) are called Hirota bilinear equations. 

Concluding this subsection let us also recall that the space of (generalized) KP Lax operators 
(|2.2| ) is endowed with bi-Hamiltonian Poisson bracket structures (another expression of its inte- 
grability) which result from the two compatible Hamiltonian structures on the algebra of pseudo- 
differential operators ^DO [49|. The latter are given by: 

{ (L\ X) , (L\ Y) }, = - (L | [X, Y]) (2.20) 
{ (L\ X) , (L\ Y) } 2 = Tt a ((LX) + LY - {XL) + YL) 

+ - fdxRes([L, xfjd^Resf^L, Yfj (2.21) 

where the following notations are used. < -|- > denotes the standard bilinear pairing in ^T>0 via 
the Adler trace (L\ X) = Tr a(LX) with Ti^X = jKesX. Here X, Y are arbitrary elements of 
the algebra of pseudo-differential operators of the form X = J2k>~oo D k Xk and similarly for Y. 
The second term on the r.h.s. of ( 2.21| ) is a Dirac bracket term originating from the second-class 
Hamiltonian constraint u r — l = on L ( p.2[ ). 

In terms of the Lax coefficient functions v r -2, ■ ■ ■ , vq, ui, v,2, ■ ■ ■ , the first Poisson bracket struc- 
ture ( |2.20| ) takes the form of an infinite-dimensional Lie algebra which is a direct sum of two subal- 
gebras spanned by {vj} and {ui}, respectively. The latter is called Wi +00 -algebra [50]. Its Cartan 
subalgebra contains the infinite set of (Poisson-)commuting KP integrals of motion Hi_\ = jTr^L" 
whose densities are expressed in terms of the r- function ( p. 10 ) as: 



dti 



ResL' 



(2.22) 



In turn, the second Poisson bracket structure ( p.21 ) spans a nonlinear (quadratic) algebra called 
Woo(r) |5l| , which is an infinite-dimensional generalization of Zamolodchikov's Wn conformal 
algebras 
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2.2 Constrained KP Hierarchy 

Definition 2.1 The function $ is called ( adjoint) eiqenfunction of the Lax operator L satis- 
fying Sato 's flow equation if its flows are given by the expressions: 

<93> is <9* k 

— =LI($) ; (2.23) 



/or £/ie infinite many times t k . 



Of course, an eigenfunction function, which in addition also satisfies the spectral equations from 
( |2.9| ), is a BA function. Arbitrary (non-BA) eigenfunctions can be represented in terms of the BA 
function: 

= / d\<p(\)i^(t, A) (2.24) 



with an appropriate "density" 4>(X) and a contour of integration T such that the integral (|2.24[) 
exists. 

Consider now the Lax operator Q = D + J2i^o ViD~ % ~ 1 satisfying the flow equations as in ( |2.3[) 
with r = 1. Let furthermore $j and ^ (i = 1, . . . , m) be, respectively, eigenfunctions and adjoint 
eigenfunctions of Q, according to ( p. 23 ). It is easy to prove the following lemma. 

Lemma 2.1 The time evolution of the pseudo-differential operator ^.D -1 ^ is given by: 

— faD- 1 *^ = [B k , SiZr 1 *,] _ ; B k = Q k + (2.25) 



Proof. The proof is a consequence of the following technical observations based on identity 
( A.2| ) from Appendix |A] and eqs.( 2.23 ): 



(B k ^ i D-H^)_ = B k {%)D-H i = ^D- 1 ^ 



and 



□ 



dt k 



(2.26) 



(2.27) 



What we learn from the above lemma is that if we define a Lax operator L, such that its purely 
pseudo-differential part is L_ = <!>Z) -1 \I', then L_ satisfies automatically the KP-type flow equa- 
tions as in ( p.3|): dL_/dt k = [B k , L]—. 

Following |21|, |22]| let d ai be a vector field, whose action on the Lax operator Q is induced by 
the (adjoint) eigenfunctions <&i,^i through: 



d ai Q = [Q , $i-D _1 ^i] 
We find from the lemma 2.1 the following proposition [21, EE 



(2.28) 



Proposition 2.1 The vector fields ("ghost" flows) d ai (2.21) commute with the isospectral flows 
of the Lax operator Q: 

[d at ,d n ]Q = n=l,2,... (2.29) 



i.e., d ai define symmetries for the corresponding Lax evolution eqs. d n Q = Q r l , Q 
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Definition 2.2 The constrained KP hierarchy (denoted as cKP r ^ m ) is obtained by identifying the 
"ghost" flow YaLi with the isospectral flow d r . 



Comparing ( 2.28| ) with the ( |2,3[ ) we find that for the Lax operator belonging to cKP rm hierarchy 
we have: 



Hence, we are led to the Lax operator L = Q r given by: 

r— 2 m 

L = L r , m = D r + J2 uiD 1 + ^D-^t 

1=0 i=l 



and subject to the Lax equation: 



8L_ 

dt k 



[(L k/r ) , , L] 



An equivalent useful form of (2.31) is: 

L = L r , m = D r + J2 u t D l + <H (D ~ k)- 1 

where 



r-2 



1=0 



i=i 



(2.30) 



(2.31) 



(2.32) 



(2.33) 



(2.34) 



It will be convenient for us to parametrize the cKP rjm hierarchy in terms of the Lax operator 
from ( 2.31 ). The functions 3>j and \Pj remain eigenfunctions of L r ^ m (|2.31| ) according to definition 
2.1 , i.e., they satisfy eqs.( |2.23|). 

We now recall results of |24|, [53f which showed that the Lax operator ( 2.3l| ) has an equivalent 
form as a ratio of two purely differential operators: 



Lemma 2.2 The cKP r]m Lax operator from \2. 3J\) can be equivalently rewritten as a ratio 



r,m — lj m+r ij m 



m, r > 1 



(2.35) 



of two purely differential operators 



(D + v m+r ) (D + u m+r _i) ■■■(D + v 1 ) 
(D + V m ) (D + V m -l) ■■■(D + v 1 ) 



with coefficients Vi,Vi subject to a constraint: 



m+r m 
3=1 1=1 



(2.36) 



(2.37) 



Proof. Let us start with the expression (2.35). It can easily be rewritten as 



r-2 



L^L- 1 = £ At I] (D + v^ 1 + A l+m+1 D l + D T 



1=1 i=l 



1=0 



(2.38) 



12 



taking into account condition ( |2.37| ). We therefore need only to concentrate on the pseudo- 
differential part of (HI) . We can rewrite it as 



m m 



-1 



ri = A ie -^' ; q m = J^ 
Let us define the quantity: 

Qi.i = (-I)*"' 



1=1 i=l 

qi = e /(«*-«*+i) 



i / t/i-i 



f\i-l+l 



, i = l,...,ro-l 



1 < Z < f < m 



Then, after using that D Qi ; i-iqi = D 1 Q\^D — Qi : i, eq.( 2.3S| ) acquires the form : 

771 777 

(L rn+r L m l ) _ = Y, rf ] I] + nD- 1 (-Qi, m -ig n 

where 



7=2 z= 



»i = n + riQi^-i i = 2,...,m 
The above process can be continued to yield the expression: 

771 

(L m+r L- 1 )_=J2$iD- 1 *i 

with 



(2.39) 
(2.40) 

(2.41) 

(2.42) 
(2.43) 

(2.44) 



7=1 



7-1 



®i = r i + E 



E 



7 — 777 7 — 777 

E E Q 



i— «i-m-a— 1 



777=1 



-l=Sj_ m _2+l S2=Sl+lsi = l 



777 Si777 



Qi— S2,i— si— lQi— si,i— 1 

, *i = (-l) m -^„ 



9m- 1 



1 < i < m 
qi (dx') m ' 1 



1 < i < m — 1 



(2.45) 
(2.46) 



We now give an opposite construction starting from the Lax operator ( 2.31 ) with pseudo- 
differential part (L)_ = J2i^=i $*-D • First, let us consider the simplest nontrivial case m = 2. 
Consider the identity: 



I 3 ^ 



for arbitrary functions (f>,ip,x> where in deriving ( j2.47| ) we have used the simple Wronskian identity 
(|A.9|) from Appendix. Here and in what follows we shall use the notations: 



W k = W k [if> u ...M = det d % ~^ 

l,J=l,...,K 

where Wk denotes the Wronskian determinant of {tpi, • • • , iftk}- 



, W = l 



(2.48) 



13 



We obtain from ( p. 47 ): 



A 



(2) 



,(2) 



-$i 



B?>=^ln W[ * 2 '* l] 



This construction can be extended by induction to yield for general m: 



(2.49) 
(2.50) 
(2.51) 



L = L r . m = L + + J2<5>iD ^ 
i=l 

= D r + 2 u*B< + £ A<f n) (D + B\ m) ) 1 (D + B^\) ^■■■{d + BM) ' (2.52) 



i=l 



with 



A 
B 



M 

(m) 



f_l^-feV$ Wm - fc+1 [^m,--.,^fc+l,^ a ] 
1 J ^ S W r m-fc[*m,...,*fc+l] 
W m _ fc+ i [^ m ,...,* fc+ i,* fc ] 



cL In 



W m _ k [* m , . . . , *fc+i] 



(2.53) 
(2.54) 



i.e., the expression ( 2.52 ) is of the same form as fl2.38| ). Comparing the latter two equations we see 
that Vi = B^ m) . a 

Let {"01, . . . , ipm} be a basis of linearly independent solutions of the m-th order differential 
equation Kip = 0, where K = (D + v m ) (D + v m -i) ■ ■ ■ (D + v{). Then one can show that: 



v k = d[ In 



W k = W h [ih,---,if>k\ 



(2.55) 



Hence by comparing ( |2.55| ) and ( 2.54j ) we obtain the following corollary: 



Corollary 2.1 The adjoint eigenf unctions *j are in the kernel of the purely differential operator 
{L m )* from 

Similarly to the Miura-like transformation fl2.53| )-( ^.54| ), we can express the Lax coefficients 
ui,A^\b^ of ( |2.52| ) in terms of the Lax coefficients Vj,vi of L r ,m in the form ( 2.35 )-( 2.36 ) as 
follows H| : 

m+r—s n s 

( d + V n s +s - V n J ( d + v n s -!+s-l ~ Vn^A X ■ ■ ■ X 

n 3 =l n s _i=l 
n3 ri2 m+r 

X (8 + V n2+2 ~ V n2 ) ^ ( 9 + V ni+1 ~ «m) ( ("» ~ U ' 



ri2=l 



ni=l 



.4 



(m) 



m+r—s i 



;=ni+l 
for s = 1, . . . , 1 — 1 
for s = r, . . . ,m + r — 1 



(2.56) 



(here it is understood that Vj = for j > m + 1 ) . 
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We conclude this section with some remarks. 



Remark. Let us point out that there is an alternative way (w.r.t. definition 22) to introduce the 
cKP r m hierarchies. Namely, one can view cKP rm in either of their equivalent forms ( j2.31| ), ( |2.35| )- 
( 2.36j ) or (2.52), as consistent Hamiltonian reductions | 25| , 28, ^] of the original full (unconstrained) 
KP hierarchy fl2.2j ) w.r.t. the bi-Hamiltonian Poisson structures ( |2.20[ )-( 2.21 ). Moreover, L = L r>m 
obeys the same fundamental Poisson bracket algebras ( 2.2C| )-( |2.21| ) as the original unconstrained 

lQ. 

Remark. The first Hamiltonian structure ( [2.20 ) for L = L r ^ m ( |2.33| ) is equivalent to the following 
free field fundamental Poisson brackets: 



{ <&i(x) , ^j{x') } x = -SijS(x - x') , rest = 
or, equivalently, in terms of the fields aj, hi fl2.34j ) : 

{ aj(x) , bj(x') } l = —5ijd x 5(x — x') , rest = 
Expanding ( 2.33j ) in a standard pseudo-differential operator series: 

CO CO 

L r , m = L + + £ U k ($, *) D~ k = L+ + ]T U k (a, b) D~ k 



k=l 



k=l 
m 



U k (*,V) = J2M-d x ) k ~^ , U k (a,b) =J2<- l ) k ~ lp k-i(- b *) 



(2.57) 
(2.58) 

(2.59) 
(2.60) 



i=l 



i=l 



where -Pfc(-) are the familiar Faa di Bruno polynomials, we obtain a series of free-field realizations 
of the Wi +00 algebra of the Lax coefficient fields ( [2.60 ) in terms of ([2.57 ) or a^bi ( 2.34| ), 

respectively. 

Remark. The second Hamiltonian structure (2.21) for L = L r>m ( |2,35[ )-( ^.36 ) is equivalent to 
the following free-field fundamental Poisson brackets: 



{ Vi(x) , Vj(x') } 2 = Uij - ^\ d x 5(x - x') , i,j 



,m + r 



{V k (x) , V t (x') }r 



{ Vi(x) , vi{x') } 2 = -d x 5(x - x ) 



hi + ~ ) d x S(x - x') 



k, I = 1, . . . , m 



(2.61) 



which, as demonstrated in refs.|27], ^9j, is precisely the Cartan subalgebra of the graded SL(m+r, m) 
Kac- Moody algebra. This latter property justifies the alternative name of the cKP rm hierarchies - 
SL(m + r, m) KP-KdV hierarchies generalizing the usual SL(r) KdV hierarchies defined through 
purely differential Lax operators: 



r-2 



L = L T , = D r + Y, u t D l = (D — v r ) ■ ■ ■ (D — vi) 



1=0 



Expanding ( 2.35| )-( |2,36j ) in a standard pseudo-differential operator series: 



r-2 



L r , m = D r + J2 <v,v)D l + ]T Uj(v,v)D- j 

1=0 3=1 

min(m— 1 J — 1) 

Uj(v,v) = £ {-iy- l - S A^lAv,v)P\ S _\ l } s (Vr, 



(2.62) 

(2.63) 
(2.64) 
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where ui(v,v), Ai m '(v,v) are given by ( |2.56| ) and 



P« (&,..., Ct> s ) = J2 ( d + ■■■( 9 + 

iiH Hs=n 



(2.65) 



are the multiple Faa di Bruno polynomials, we obtain a series of free-field realizations of the 
nonlinear W^) algebra of the Lax coefficient fields ui(v,v),Uj(v,v) ( |2.56j ),( pl34| ) in terms of Vi,vi 



(2.61) 



3 Darboux-Backlund Solutions of Constrained KP Hierarchies 

3.1 DB Transformation and Eigenfunctions of the Lax Operator 

We begin with reviewing the basic properties of DB transformations for the general (unconstrained) 
KP hierarchy stressing its basic property of preserving the form of the pertinent Lax evolution 
equation. Related material can be found in, e.g., refs.[54, ^jj ]. Next, we discuss in more detail the 
DB transformations and explicit solutions for the constrained cKP rjm hierarchies. 

To describe a result of acting with the DB transformation in the setting of pseudo-differential 
operators we need first to make a simple technical observation. 

Lemma 3.1 For an arbitrary pseudo- differential operator A we have the following identity ]5^/ : 
'xDx^AxD^x' 1 ) = xDx^A+xD-'x- 1 - X d x (V^+Cx)) ^V* (3.1) 



Let L satisfy the generalized KP Lax equation d k L 



. Consider a "gauge" transformation: 

(3.2) 



L-> L = TLT~ 

with T being a pseudo-differential operator. Then the transformed Lax operator will satisfy: 



duL 



TL^T- 1 + {d k T)T~ l , L 



(3.3) 



Definition 3.1 The DB transformation is defined by ( \3. with the following special form ofT : 



-i 



T = $D<S> 

i.e., with <3? being an eigenfunction of L. 



One easily verifies that for A = L~ and x = $ eq. (|3.1| ) becomes: 



(3.4) 



TL~T 



Correspondingly, (|3,3j) takes the form: 



-i 



d k L 



TL r jT~ + (d k T)T- 



Lv,L 



(3.5) 



(3.6) 



and we have, therefore, established: 



16 



Proposition 3.1 The DB transformation \3.Q) with an eigenfunction $ (3.4 ) preserves the form of 
the Lax equation (2^5), i.e., the DB transformed Lax operator satisfies the same evolution equation 
(3.6) as the original Lax operator. 

Using simple identities valid for any pseudo-differential operator A and arbitrary function / : 
Res (/A/ -1 ) = Res A , Res (dAD^ = Res A + d x A {0) (3.7) 
where the subscript (0) indicates taking the zeroth order term, we obtain: 

ResZ- = Res (qD^L^D- 1 ®- 1 ^) = ResL^ +a r -^-ln$ (3.8) 

Combining ([O]) with eq.( |2.22| ) for the r- function, we get therefore 54]: 



Proposition 3.2 Under DB transformation $3. 1 ) with an eigenfunction <3? (3-4) the t -function 
associated with the Lax operator L transforms according to: r — » r = $r. 

It is also useful to write down the DB transformation in terms of the dressing operator and the 
BA function as well: 

W -> W = ($£>$ _1 ) WD' 1 (3.9) 
W, X) - m X) = \-i*(t)d ( 2 ^) (3-10) 

Inserting (|3.10| ) and f(t) = ®(t)r(t) (from Prop.U) in the expression fl2,lC| ) for the DB-transformed 
BA function: 

m\) = fit Z^ 1]) e^ (3.11) 

T[t) 

we deduce the following identity relating the r-function with an arbitrary eigenfunction of L : 

*(' -If'D - ! + A-.flm.t,) _ A -i flh T C - If!) (3.12) 

r(t) v y 

Let us illustrate possible applications of ( |3.12| ) by the following example. Taking into consideration 
that: 

A»"-'f» = £ ^-'fl'fj" (3.13) 



and expanding in A we find from (|3.12j ) : 

rit-iX- 1 ]) _ P2 (-[dmt) 



r(t) A*(t) 



+ 0(A^ 2 ) (3.14) 



This implies that: 



0nr) »_B(-iai)*_i(a-ft)« ( 3.i5) 

v y $ 2 $ 



1 ((^ 2/r ) + " 5?) $ _ i 



2 $ r 



-■u r _ 2 (3.16) 
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Hence by only using the basic fact that $ was an arbitrary eigenfunction we obtained from (3.12) 
a very general result (lnr)" = ResL 1//r . 

Analogously to fl3.12| ) one also obtains an identity relating the r-function with an arbitrary 
adjoint eigenfunction <I> : 

* {t + [A " ]) - 1 - A-^lntfCt) = -A-^ln (3.17) 



*(*) 



T{t) 



The next step is to consider DB transformations within cKP rm hierarchy, i.e., such transfor- 
mations ( |3.2| ),(|3~i|) which preserve the constrained cKP ri?n form ( 2.3 1| ) of the corresponding Lax 
operator. Under an arbitrary DB transformation L = (x^X 1 ) L (xD~ X )i where x is an 
eigenfunction of L ( p. 31 ), the transformed Lax operator reads: 



L 



La 



X 



i=l 

L+X 



D~ 



>i 



o 



i=l 

m 



(3.18) 
(3.19) 
(3.20) 



X 



i=l 



x£>x _1 £) (x) 



(3.21) 
(3.22) 



Let us recall, as already stressed in the introduction, that all functions appearing in ( |3.19| )-( |3.22| ) 
are (adjoint) eigenfunctions of L ( [2.31 ) or L ( |3.19| ), but they are not (adjoint) BA (wave) functions. 
In other words, they satisfy: 



d_ 

dt k 



f = Ll(f) f = X ^i 



9 T 

dt k 



(3.23) 



and similarly for / = <3>o, &i , \Pi with L — > L. 

In case when x coincides with one of the original eigenfunctions of L, e.g., x = $i> it follows that 
$1 = and the DB transformation (3.18) preserves the form (2.31) of the Lax operators involved, 
i.e., it becomes an outo-Backlund transformation. Applying successive DB transformations 



L (k) _ T (k-i) L (k-i) ^ r (fc-i)^ 1 



yields: 



(fc+1) 



$ (*+l) = ^ T{ k) L {k)^ ^(fc) 



KT 1 



fc = 0,1, 



, m 



(fc+i) 



T (fe)\ 



1 



(*) 



-i 



l /\t,(*)*(*) 



alt 1 *r*r 



(3.24) 

(3.25) 

(3.26) 
(3.27) 
(3.28) 
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Using the first equality from fT25p , i.e., L^+^tM = T^L^ , one can rewrite ( |3.26| ) in the 
form: 

(3.29) 
(3.30) 



whereas: 



$ (fc) = r (fc-i) T (fc-2)... r (o) ((^ (0) ) ($i 0) ) 

, i = 2, . . . , m 



= T (fc-l) r (fe-2) . . . T (0)^(0)^ 
Accordingly, for the BA functions we have from ( |3.10| ) : 

V> (fc) (i, A) = A -1 ^*- 1 )^^- 1 ^*, A)) = \-k T {k-l) T (k~2) . . . T (0)^(0)^ A )) ( 3 _ 31 ) 

which together with ( |3.29 ) and ( |2.24 ) implies the following "spectral" representations: 



-( k \t) 



= J dXX k cf>f\x)^ k \t,X) , i = 2,..., 



m 



(3.32) 
(3.33) 



where A) was defined in 

Finally, for the coefficient of the next-to-leading differential term in ( p. 31 ) u r ._2 = r Res L~ = 
r dollar (recall eq.( 2.22j )) we easily obtain from fl3,19| ) (with x = hs A;-step DB-transformed 
expression: 

T (fc) 



( fc ) _„(°) 

-2 "r-2 



(0) 



which conforms with the result of multiple application of Prop. 3.2 



_(*) 

i(oy 



(3.34) 



(3.35) 



Let us consider momentarily the special case of cKP r i hierarchies (i.e., m = 1). Combining 
eqs. dj2|),0) (for m = 1) and flpg ), we deduce the following: 



Proposition 3.3 ^4ny Darboux-Backlund orbit of cKP r i hierarchy (2.31) (for any r > 2 and 
m = l,) defines a structure of two-dimensional Toda lattice model j^j (here <E> = Q\) : 



x -^-ln*( n > = — ^ 

dt r 



n = 0,1,2, 



-i 



5 



d x — lnrW 



r (n+l) T (n-l) 



(tW)' 



, 71=1,2,... 



(3.36) 



(3.37) 



Indeed, introducing variables such that 
-1,-2,... , the recurrence relation (3.36) for 



$( n ) = expjV'n+i — ip n } and ^ = for n = 
the cKP r i eigenfunctions turns into the familiar 



two-dimensional Toda lattice equations of motion, whereas the recurrence relation eq.( |3.37 ) for the 
corresponding cKP r i r-function r' n '(t) coincides with the equation for a (partial) two-dimensional 
Toda lattice r-function. 

For r = 1 case (recall t\ = x) eqs. (|3.36| ),( |3.37|) degenerate into one-dimensional Toda lattice 
structure (see eqs.(^8^),( ^^ ),(^^) below). 
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Returning to the general cKP rm case, we can represent the fc-step DB transformation (3.29) 
( 3.35| ) in terms of Wronskian determinants involving the coefficient functions of the "initial" cKP ri 
Lax operator: 



r— 2 m 

L (0) = D r + £ uf ] D l + Y, gWd- 1 *® 



i=l 



only. Defining: 



= 1,2,... 



(3.38) 
(3.39) 



we arrive at the following general result: 

Proposition 3.4 The k-step DB-transformed eigenf unctions and the tau-function ft3.2!\ ) — (3.31 ) 
of the cKP rm hierarchy for arbitrary initial ( 3.38 ) are given by: 



(k) 



(k) 



, J 



, m 



-(*) 



W k [^\ X {1 \...,X (k - l) ]r 



(fc-i)w(o) 



(3.40) 

(3.41) 
(3.42) 



where t^°\t^ are the r -functions of L^°',L^ k \ respectively, and x® *s given by ( 3.3Q) . 

The proof of ( |3.40| )- (|3.42| ) is accomplished via induction w.r.t. k by using the important composition 
formula for Wronskians ( |A.9| ) ( [A. 1 1|) from Appendix [A]. 

Now, comparing expressions (|3.32| ) and ( g!4C| ) together with (|1|), we get the following recur- 
rence relation for the T-functions of the DB-orbit on cKP rm : 



.(fe+i) 



(*) 



d\X k ^ ( f>f ) (X)e^T^(t-[X- 1 ]) 

' dX A fc ( 1+ ^ 0) (A) : e~ § ^ : (t) (3.43) 

where in the second eq. fl3.43| ) we have employed the well-known vertex operator |6|, |30], ^0| with: 

'«> -VA + V-,-- (3.44) 



where : . . . : indicates standard normal ordering, meaning J^- "modes" to the right of t\ "modes" . 
Note, that eq.( 3.43| ) is nothing but a special case of Sato Backhand transformation for the con- 
strained cKP ri?n hierarchies. Using the Wick-theorem identity: 



: e 



-0(A fc ) 



: e 



-0(Ao) 



=: e 



(3.45) 



z=o 



i>3 



one can solve the recurrence relation ( |3.43j ) expressing r( fc+1 ) (t) directly in terms of the "initial" 
r-function t<® of L<°) (KM : 



.(fc+i) 



(<) 



k k 

IJ dX j W k+1 [f (t, X ),d r f (t, AO, ... , d k f (t, A fc )] rW(i - E^ 1 ]) ( 3 - 46 ) 
j=o i=o 
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tiv [u dA i (n & - (n ( A r - A rt) n m*,a s ) ^(t - bvd ( 3 - 47 ) 

" + " ^ j=0 i>j i>j s=0 j=0 

where we used the simple antisymmetrization identity: 

/ N \ 



(3.48) 



antisymm {Wj m} 



\S-1 



V*=l 



1 

iV! 



n - us) 



(3.49) 



Now, comparing expressions ( 3.46 ) and Q3.42|) for the Wronskian r- functions we find that they 
agree provided the following identity involving r(°) is satisfied: 



n - a j) 

i>j 



^(t-EtoCA- 1 ]) 



r(°)(t) 
^^-[Ao 1 ]) 



(3.50) 



-30,3X,—,Jk 



(0 + a*: 



r(°)(i) 



(9 + A 



■ (0) (*-[A^]) 



-(°)(t) 



It is easy to see that eqs.( |3.50| ) coincide with the well-known differential Fay identity (for k = 1) 
and its higher derivative corollaries |4C]. 

One can generalize the successive DB transformations ( |3.29D on general cKP rim Lax operators 
( 2.31 ) as follows. Within each subset of m successive steps we can perform the DB transformations 
w.r.t. the m different eigenfunctions of ( p. 31 ) unlike ( p.29 ) where all DB transformations are given 
by the first $i eigenfunction solely. Repeated use of eqs, ([A,9|) -( A.10| ) from Appendix and employing 
short-hand notations: 



X 



1, . . . , m 



(3.51) 



yields the following generalization of ( 3,40|) - (|3,42| ) 



Proposition 3.5 The most general discrete Darboux-Backlund orbit on the space of the cKP r>m Lax 
operators, starting from an arbitrary initial (3.3t), consists of the following elements: 



r—2 m 
L (km+l) = jy + u ^ km+l) D j + y] $( fcm+i )£)- 1 ^( fcm+ ') 

j=o J i=l 



(3.52) 



with k arbitrary non-negative, 1 < I < m, and where: 



W 



(km+l) 



-,(km+l-l) 



-I 



T ,(fcm) T ,(fc m -l) nn((Jfc-l)m) rp(m-X) rp(P) (k±) 

. . . j. 1 i m . . . i 1 . . . x m . . . i l Xi 

a(0) (!) (!) C*- 1 ) ( fc ) ( fc ) ( fc ±) 

$\ J ,...,$W,Xl ,---,Xm,---,Xl ,---,Xm ',Xl ,---,Xl ,Xi 



w 



*(°) *(°) V (1) Y (1) Y (fc_1) V (fc_1) Y W Y W 

, . . . , , Xi ; ■ ■ ■ > Xm > ■ • • , Xi Xm > Ai > • • • > a; 



(3.53) 



x ( fc +) = x ( fc +!) fo r 1 < i < / 



' = for / + 1 < i < m 



The corresponding r functions read: 
T (km+l) 



T 



(0) 



(fcm+Z-1) 
I 



§(km) $ (fem-l) *((*-l)»n $(m-l) $ 
• " " 1 m • " " ^ 1 ■•' m ...a; 



(0) 



W 



(0) 

i ) 



$(0) y (l) y (l) 

T m lAl ! • • • j Am > ' ' ' i Al 7 



v (*-l) y W 
5 Am i Al ) 



<Xi 



(fe) 



(3.54) 
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3.2 Truncated KP Hierarchies 



An important particular class of cKP rm DB orbits ( p.52 ) are those with a purely differential initial 
Lax operator ( 3.38] ) : 



r-2 



L(°) = (l(°)) =D r + J2^f } D j , i.e. * 



(0) 



j=0 



(3.55) 



(k) 



Such cKP r>m DB orbits are characterized by the fact that the m adjoint eigenfunctions ^ = ^ 

(k) 

are not independent from the m eigenfunctions <£j = &\ since both are parametrized in terms of 
m independent initial eigenfunctions <E>^ only. 

Formulas (3.40)-(3.42), ( 3.53] )-( f3.54 ) simplify significantly in this case since now we have: 



x w s ( L( o)y^ = w (o) fi = li 

leading to (e.g., for (g!4C|)-(g!4^)) : 

W k+l [^\d r ^\...,d^} 

W k+l [^\d^\...,d^\*f) 



, m 



(*) 
l 

(*) 

3 

T (k) 

T 



(0) 



j = 2, . . . , m 



(3.56) 

(3.57) 

(3.58) 
(3.59) 



A special feature of the subclass of cKP rjm DB orbits with a "free" initial = D r is that 
their dressing operators are truncated (having only finite number of terms in the pseudo-differential 
expansion, cf. third ref. in pl| ) : 



yy{km+l) 



T, 



(km+l-1) r r (km)r r (km-l) rpdk-l)™,) rpfm-l) rp(0) 7-)- 

l ' ' ' J 1 1 rn ' ' ' 1 l " ' 1 m ■ ■ ■ ± 1 LJ 

km+l 

wf m+l) D-i 

j=0 



■km— I 



(3.60) 



where notations fl3.51|) were used. Recalling fl2.14j ), eq.( 3.60| ) implies: 



{km+l) = Pl (-[d])T^)(t) 
i T (fcm+i)(t) 



Pj{-[d})r {km+l) (t) = , j > km + 1 + 1 (3.61) 



On the other hand, one can consider a priori truncated W dressing operators outside the context 
of Darboux-Backlund transformations. 

Definition 3.2 KP hierarchies defined through Lax operators built out of m-truncated dressing 
operators: 



W = W (m) = Y j w i {t)D- 



L = L< m > 



i=0 



W {m) D (V (m) ) 1 



_d_ 

dt„ 



w {m) = - (w^ m) D(w {m) y l ^j 



(3.62) 
(3.63) 



are called m-truncated KP hierarchies. 
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As in flirftiD ab ove, the r-function of the m-truncated KP hierarchy obeys the constraints: 

Pj(-[d])T^{t) = , j > m + 1 (3.64) 

Let us set m = 1 and recall the Hopf-Cole transformation, which describes the dressing and Lax 
operators by a function </> such that w% = —d x ln(j) (see eq.( 3.82j ) below). Hence <fi is a solution of 
W^90 = 0. Similarly, for arbitrary m we generalize the Hopf-Cole transformation, and consider 
a set of m functions being solutions to the system of differential equations: 



W im) d m (/) k = k = 1, 



(3.65) 



It is well-known, that while solutions of the general KP hierarchy form the universal Grassmann 
manifold UGM, solutions of ( |3.65D defining the m-truncated KP hierarchy form the Grassmann 
manifold GM(m, oo) = Mat(oo x m)/GL(m; (TJ) where Mat(oo x m) denotes oo x m matrices of 
rank m |56|, |57| . 

In terms of the solutions of the m-th order differential equation ( |3.65| ) the Wilson-Sato equations 
( 3.63| ) take the simple form: 

d n 4>i = d^cpi i = l,...,m (3.66) 



Note that fa satisfying ( 3.66| ) can be regarded as a set of m independent eigenfunctions of the "free" 
Lax operator £A°) = D, and their explicit form reads: 



((*) 



(3.67) 



i>i 



with arbitrary "densities" <^°^(A). Moreover, it is assumed that in general the set {(pi, 
nondegenerate in a sense that no one fa is a derivative of another 
In different words we have the following Lemma: 



i} is 



3 ' 



Lemma 3.2 Eqs. ( 1T6(\) are equivalent to the Wilson-Sato equations ( 3. 6$) for the m-truncated 



dressing operator (3.6i). 



Proof. Define for convenience T = WD m = D m + . . . , which is a purely differential operator of 
order m. The Lax operator can then be rewritten as L = TDJ 7 ^ 1 (here we will be suppressing for 
brevity the superscripts ( m ) on W and L). Assume that ( 3.66[ ) holds. It follows: 



= d n {F<j> k ) = {d n T)<j) k + F{d r x 



(d n T + L n _F) 



(3.68) 



for some arbitrary integer n. In (3. 



use was made of the obvious identity L n = TD n T~ x or 
U\_ T + L'^J 7 and the fact that L^Tcpk = 0. Expression {d n !F + T) on the right hand 
side of ( |3. 68 ) is a purely differential operator of order smaller than m. Since by assumption all m 
functions <f>k are independent we conclude, therefore, that d n T = —U!lT, which is equivalent to 
the Wilson-Sato equation ( |2.5D for the dressing operator W. 

Let us now prove the inverse statement. From assumption d n T = L+J 7 — TD n and ( 3.65| ) we 
find that J- {d n (f) k — d™(j) k ) = 0. Since (fi k span the m-dimensional space of solutions of T operator, 
we must have d n (j> k — d™(t>k = J2j CfcjVj'j with c$ j,k = 1, . . . , m being x-independent constants. 



Hence, V n 0fc = d^cp k for the "covariant" derivative (V 
then [ V„ , V; ] <\> k = leading to c[j' 
4>k = 4>l ( ex P c )ifc- It follows that d n (f>k = d™4>k- D 



» 



jk = ^nOjk - c kj . Compatibility gives 
c~~ k d n Ckj for some x-independent matrix cy. Define now 
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Let us return to (3.65) which factorizes as follows: 



Ffa = (D m + Wl D m ~' + ■ ■ ■ + w m U fc = (D + v m ) {D + Vm -i) ■••(D + v 1 ) 







(3.69) 



There is a relation between coefficients V{ of ( p. 69 ) and solutions 4>k of the same form as ( |2.55| ) with 
<j>k replacing ip k : 

Wi-x\6 x , . . . , 



9 hi- 



i 9i 



) • • • > vtj 



1 



(3.70) 



and, therefore: 



•Ti 



--P — = D + Vi 



Wi-i Wi 

L (m) = W {m) D ( W {m)yl = J^ D jr-l =T m ---T 1 D Tf 1 • • • T" 1 



(3.71) 



Comparing (3.71) with ( p. 25 ) we identify the m-truncated KP hierarchy ( 3.62| )-( |3.63| ) as the m-th 
member of the DB orbit of cKPi m with initial ZA°) = D. Accordingly, we have the following: 

Proposition 3.6 The r-function and the BA function of the m-truncated KP hierarchy take the 
following explicit form: 



_(m) 



p s {-[9])W m [tpi, ...,4> r, 



w, 



m \a>i, . . . , <p m \ 



det 



d 



i-l. 



0, 



i,j=l,...,m 

for 1 < s < m 
for s > m + 1 

wt- is \ 



s=l ii<...<i s 
, . . . , d(p h d(p is ,...,(pm[ 



W„ 



(3.72) 
(3.73) 

(3.74) 
(3.75) 



where the Wronskian (3.75) is obtained from W m = W m [4>i , . . . , 4> m ] by putting an additional deriva- 
tive d on rows with numbers i\, . . . ,i s , respectively. In particular, we have: 



Pk(-[d])4>i = for k > 1 , i = l,...,m 



(3.76) 



which is identical with (3.6t). 
Accordingly, eq.( [3.62 ) now reads: 



W (m) _ 

3=0 



it 



[ D -3 



(3.77) 



m [<pi, . . . , (p m \ 



Proof. Due to the above identification of m-truncated KP hierarchies, eq.( |3.72| ) follows straight- 
forwardly from (|3.59|) with r = 1. Furthermore, we have: 



W m [(/>i(t - [A -1 ]), ...,4> m (t— [A -1 ])] = ^ A _s p s ( - [<9] ) W m [(f>i (t),...,(f> m (t )] (3.78) 

(3.79) 



s=0 
m 



£(-l) s A~ s W t~ is) 

s=0 ii<...<i a 
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where use is made, respectively, of the definition of Schur polynomials ( |2.13| ) and of the simple 
identity (cf. (|j37|)) : 



(3.80) 



which follows from e xp f^/l^ l j = (1 — A 4 / A) • Then, eg. ( 3. 73; ) directly follows upon com- 

paring of (|3.78|) with ( 3.79 ). Accordingly, eq.( 3.74 ) is easily verified upon substituting fl3.72|) in the 
expression for the BA function: 



*W(t A) = T(m)(t - [A ' 1]) e^) 
V r M(t) e 



(3.81) 



and using again ( 3.79| ). □ 

3.3 Generalized Burgers-Hopf Hierarchies 

For the simplest case m = 1 relation ( |3.69| ) takes the form of a classical Hopf-Cole transformation: 
(l + wiD" 1 ) D(<f>) = -» (a + ioi)^» = -> Wi = -S a! ln^ (3.82) 
Correspondingly, T = W^D = D — d x (ln</>) = cfyDcj) -1 and, thus, we are lead to the Lax operator: 



D + 



4>(hx4>)' 



D 



(3.83) 



One finds from Lem ma |3.2j that the Lax equatio n (p.3[ ) for LW is equivalent to d n (j) = d™(f>, which 
upon comparison of (|3.830 with (]3.251) identifies fl3.83|) as k = 1 member of cKP^i DB orbit. 

Equations d n (j> = d™(f) can be rewritten in terms of the coefficient w = —w\ of the dressing 
operator = 1 — wD~ l as: 



d n w = d x P n (w) 
P n+1 {w) = (d + w)P n {w) n = 0,1,2, 



PoH = i 



(3.84) 
(3.85) 



where P n (w) are Faa di Bruno polynomials fully determined by the recurrence relation in ( J3.85Q . 
The system of nonlinear differential equations ( 3.84 ) , which is an alternative form of the equations 
of the k = 1 member of cKPii DB orbit, is called a Burgers-Hopf hierarchy. 

Let us now consider the fully degenerated m-truncated KP hierarchy, i.e., for which: 



,3 = 1, 



, m 



(3.86) 



Substituting djgg ) in (|]7|) and comparing with ( |3.59| ) allows us to identify it with the m-th 
member (m > 1) of the cKPii DB orbit starting from = D. In analogy with the above 
identification for the (ordinary) Burgers-Hopf hierarchy, we shall introduce the following: 



Definition 3.3 We call the fully degenerated m-truncated KP hierarchies generalized Burgers-Hopf 
hierarchies. 
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Alternatively, the generalized Burgers-Hopf hierarchy is described as a Darboux-Backlund orbit of 
the Burgers-Hopf hierarchy defined in ( 3,83 ). Below we demonstrate that the set of generalized 
Burgers-Hopf hierarchies is also equivalent to the semi-infinite one-dimensional Toda chain system. 

The Lax structure obtained in the process of successive DB transformations applied to Burgers- 
Hopf hierarchy ( 3.83 ) takes the following form of recursive relations: 

L {k+i) = ($(k) D $(k)- ly ) D-^r 1 ) = D + $( fc+1 ) D~ l ^ k+1 > (3.87) 

= $(fc) An^WV' + ($(k)\ 2 y(k) ; ^(fe+i) = ^(fe)V 1 ( 3 88 ) 

with $(°) = cf). 

Proposition 3.6 specializes in the case of generalized Burgers-Hopf hierarchies (3.87) to the 
following: 

Proposition 3.7 The r-function and the BA function of generalized Burgers-Hopf hierarchies 
@J$) read: 



J 3=1 i>3 «=1 

Ps{-[d\)W k [^d^...,d k - 1 <P] = 

_ f {-l) s W k [(t>,d(f ) ,...,d k - s - 1 (f>,d k - s+1 4>,...,d k 4>} , for 1 < s < k 

for s > k + 1 

WM, dcf),..., d k - s - 1 ^, d k - s+1 (p, ...,d k ^ 







^ fc )(t,A) = e^){l + 5:(-A) 
$(*)(t) 



w k+1 [(i>,d<!>, ■■■,&</>] 



(3.89) 

(3.90) 
(3.91) 
(3.92) 



The second equality in ( 3.89| ) is a special case of eq.( 3.47| ) with r = 1, r' ) = const . 

Using the second equality in ( 3.79| ) in the fully degenerated case of Q) we also find further 
useful identities involving Wronskians of the type appearing in (3.91) : 

J d\X 2k ^^°\X)e^W k [ ( l ) (t- [A" 1 ]),^- [A- 1 ]),...,^" 1 ^- [A" 1 ])] 



ifc— s— 1 J, ofc— s+1 



f Wjfc+x^a^,...,^] , for j = 

I , for I < j < k 



,9 



k-l 



(3.93) 



In particular, identity ( |3.93| ) with j = provides explicit solution for the r-function recurrence 
relation ( 3.43| ) (with r = 1). 

In order to exhibit the semi- infinite Toda chain structure in fl3.87D - fl3.88D , let us introduce the 
following change of variables: 



k = 0,1,... 



(3.94) 
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Substituting ( 3.94| ) in ( |3.88| ), we find that \jj n satisfy the standard form of the one-dimensional Toda 
lattice equations: 



^Ipn+l+lpn-l—Zlpn 



(3.95) 



with if) n = for n < 0. Furthermore, comparing flpj) with (g^) implies that eq.(ggD has 
solution in the form of a Wronskian: 



ip n = In W n [4>,d(f>, ... ,d 



,n—l . 



with 



,01 



(3.96) 



Let us return to the constraint equation ( |3.76| ) in the context of generalized Burgers-Hopf hier- 
archies (fully degenerate m-truncated KP hierarchies) : 

p k (-[d])<f> = for k > 1 (3.97) 

which is identical with ( 3.66| ). There exists a solution of this condition which can be written as: 



K 



EOO . I 



(3.98) 



r=l 



with arbitrary constants c r , X r , r = 1,...,K. Alternatively, the solution ( 3.98 ) can be directly 
obtained from the "spectral" representation ( 3.67| ) by making the particular choice |)8| : 0^°^(A) = 
J2^=i c r 5{\ — A r ). One can easily prove that: 



Lemma 3.3 



K 2 K 

w K [<t>, d4>, ... , d K - vi = n c r ii( A i - ^wkKu c 2 , . . . , (k) = (n (^» - A i)) n ^crw (3.99) 



r=l i>j 



r=l 



Eq.(3.99) establishes relation between, on one hand, the DB construction of the r-function and, on 
the other hand, the Wronskian type multi-soliton solutions of the nonlinear evolution equations of 
the (constrained) KP hierarchy common in the literature. Hence, the solution space for the semi- 
infinite one-dimensional Toda chain formulated as a generalized Burgers-Hopf system contains the 
well-known multi-soliton solutions of the constrained KP systems. 



4 Additional Symmetries of Constrained KP Hierarchies 
4.1 Background on Additional Symmetries 

Let L be again as in eq.(|2.2|) a pseudo-differential Lax operator of the full generalized KP hierarchy: 



r-2 



L = D r + J2 VjD j + Y, u i D ~ { 
j=0 



i>l 



d_ 

dti 



L 



LI, L 



(4.1) 



(recall x = t\) and let M be a pseudo-differential operator "canonically conjugated" to L such that: 





L , M 



t 



dti 



M 



L r + , M 



(4.2) 



27 



Within the Sato- Wilson dressing operator formalism, the M-operator can be expressed in terms 
of dressing of the "bare" operator: 



1>1 l>l 1=1 

conjugated to the "bare" Lax operator L^ ' = D r . The dressing gives 

M = WM^W- 1 = WX (r) W- 1 + l ~^t r+l Lr = l ^^t r+l L\_ + M_ 

i>i r i>o r 

M_ = WX^W- 1 -U-E — • W- 1 



(4.3) 



(4.4) 
(4.5) 



where in (|4.5|) we used eqs. fl2.5|) . Note that Xr r -\ is a pseudo-differential operator satisfying 



On BA functions ( |2.9| ) and eigenfunctions (2. 23), (2. 24) of L the action of M is as follows: 

f \ l - T d 



Mif;(t, A) 
M$(t) 



+ a r (X) i>(t,\) 



dX 



d_ 

dX 



+ rX r ~ 1 a r {X) 



X 



l-r 



-<t>(\)U{t,X) 



(4.6) 
(4.7) 



where a T (X) is a function of A onlyQ. 

The so called additional (non-isospectral) symmetries 30 1 are defined as vector fields on the 
space of KP Lax operators (4.1) or, alternatively, on the dressing operator through their flows as 
follows: 



dh._n.L 



(M n L k ^j ,L] = [(M n L fe ) , L] +nM n - l L k ; d Kn W = -(M n L k ^j W (4. 



The additional symmetry flows commute with the usual KP hierarchy isospectral flows given in 
(|2.3| ). But they do not commute among themselves, instead they form a Wi +00 algebra (see 
e.g. [p0[| ). One finds that the Lie algebra of operators <9fc jn is isomorphic to the Lie algebra 
generated by — z n (d/dz) k . Especially for n = 1 this becomes an isomorphism to the Virasoro 
algebra d k l ~ -£ fc -i, witl1 [Ai , C k ] = (n- k)C k . 



4.2 Modification of Additional-Symmetry Flows for cKP r m Hierarchies 

The main task of this section is the formulation of additional non-isospectral Virasoro symmetry 
structure for the cKP rjm hierarchies, or in different words, we will redefine the standard additional 
symmetry flows (|4.8| ) for the full KP hierarchy in such a way that they will preserve the subspace 
of cKP rjm hierarchies ( |2.31| ). To be more precise, we first show that only the si (2) subalgebra 
(containing Galilean and scaling symmetries) of the standard Virasoro additional symmetry algebra, 
with flows given by dk,i for k = 0,1,2, preserves the constrained form of cKP r m Lax operators. 
Next, we show how the rest of Virasoro additional symmetry (for the Virasoro generators C k , k > 
2) acting on cKP r m is recovered through an appropriate modification of its generators. This is 



4 The appearance of a r (X) can be traced back to the ambiguity in the definition of the dressing operator (2.4): 
W — ► WWo where Wo = 1 + X<>i c ^~ l with constant coefficients d. 
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accomplished by adding a new structure consisting of "ghost" flows ( 2.28| ) related to the plethora 
of (adjoint) eigenfunctions characteristic for the cKP r]m hierarchies. This is natural since the 
"ghost" symmetry flows commute with the ordinary isospectral flows (prop. |2.1| ). 

Let us begin by applying the additional-symmetry flows (4.8) on L Q2.3ip . For n = 1 we get: 



(d k ,iL)_ 



ML' 



L 



L 



(4.9) 



Using the simple identities ( 2.25| ) and ( A.4 ) from Appendix [A] for the Lax operator ( [2.31 ), we are 
able to rewrite (fO|) as: 



(d k ,iL)_ 



8=1 

m k—1 



1=1 



i=l j=0 



Here 



(4.10) 



(4.11) 



l^^M^ + E^A -1 ^) 

(and similarly for the adjoint counterpart) denotes action of L on <!>,. Notice that L k ~j~ 1 (<fri), 
(L*Y (^j) are (adjoint) eigenfunctions of L ( |2.31| ). Hence, whereas the original L ( 2.31[) belongs to 
the class of cKP rm hierarchies, the transformed Lax operator given by <9 M L (cf. eq.flO])) belongs 
to a different class - cKP r)m (^_i) (for k > 3), since the number of eigenfunctions has increased. 

For k = 0,1,2 the flow equations ( |4.10[ ) can still be rewritten in the desired original cKP r!m 
form: 

m 



1=1 



with d T = dk,i (k = 0, 1, 2), where: 

Bo^i = (Af)+ ; 3>,i*i = - (Af)+ (*i) 

M i = (Mi;) + (* i ) + Q$ < ; ^^-(MLjl^J+^i a + /? = l 



d 2) l$i = ML 2 (*0+L($ 



a 2il ^ = - (ml 2 ) (*i)+L*{*i) 



(4.12) 

(4.13) 
(4.14) 

(4.15) 

The coefficients a, (3 on the right hand sides of ( 4.14| ) are fixed to a = f3 = \ from the requirement 
of closure of the algebra of the flows ( 4.13|) -( [L15 ). 

Since the additional flows satisfy the algebra [di t i, = — {I — k) , we have an 

isomorphism 3^,1 ~ — £k—l with the Virasoro operators and equations fl4.13| )-( ^.15| ) contain the 
si (2) subalgebra generators Cq, d\. 

However, for d T = d^i k > 3, eq.( 4.12| ) does not hold anymore due to absence of consistent 
definitions for dk,i<&i, dk,i^i generalizing ( f4,13| )-( [4.15| ) for higher k. Thus, it appears that the 
symmetry constrains behind the cKP r]m hierarchies have broken the standard Orlov-Schulman 
additional Virasoro symmetry down to the sl(2) subalgebra. 

To recover the complete Virasoro symmetry, our strategy will be to redefine the additional- 
symmetry generators. We first describe our technique for k = 3 in which case eq.( |4.9| ) contains a 
term: 



■in in m 

(l 3 ) = e ^D' 1 ( L *) 2 + E L D ' lL * (*0 + E l2 (**) D '^i 



(4.16) 



8=1 



8=1 



8=1 



29 



Note that the middle term in (4.16) is not of the form of required by equation (4.12). At this point 
we recall that for 

X = M k D- x N k (4.17) 

k=l 

with definitions (fE3lD and ( p7|) we find using identity (|A.3j ) from Appendix |A|: 

/ m 

[X, L]_ = Y / (-HM k )D- 1 N k + M k D- 1 L*(N k ))+Y / (x(^i)D~ 1 ^ i -^ l D- 1 X*^ i )) (4.18) 



k=l 



i=l 



According to Proposition |2.1| the flows generated by ( 4.18| ) will commute with the isospectral flows 
(|2.3| ) provided Mi, Ni are (adjoint) eigenfunctions, which will be the case in what follows. Consider 
now as an example: 

Y 3 = -J^^iD^L* -L^D- 1 ^,, 



i=l 



(4.19) 



Using the above formulas we find that 



o rn 

[Y 3 ,L]_ = - (L 3 )_ + - £ (^ZT 1 (L*) 2 + L 2 (<&<) D- 1 ^) 

2=1 



i=l 



(4.20) 



Hence 



(ML 3 )_+y 3 ,i 



It turns out that it is possib 



still has a form of ( 4.12 ). The interpretation of this result is evident, 
e to find additional symmetries for the cKP r m model by combining the 
original d kj \ flows and the ghost flows Q2.2S| ) associated with operators of Y 3 type. This will work 
provided that the above construction yields the Virasoro generator £2 obeying the correct algebra 
with the unbroken si (2) generators found above in ( (4.13|) ~( |4.15| ). Note that each of the two terms in 
( 4.19j ) could have been used with an appropriate factor to obtain the similar conclusion as in (|4.20j ). 
The choice of the coefficients in Y3 ( 4.1S| ), apart from being the most symmetric combination, has 
the advantage that it will lead below to the correct Virasoro algebra commutators. 

We now generalize the above manipulations to an arbitrary k. We introduce the pseudo- 
differential operators: 



4 0) 


m k—1 

= EE 




i=l j=0 




m k—1 


- EE 




i=l j=0 




m k—1 


- EE 



k > 1 



i=i j=o 



j(k-l) + 



k > 1 



(jb-2)(Jfe-l) 
6 



(L*y (*j 



(4.21) 
(4.22) 
k > 1(4.23) 



Note that x!p = Y 3 as defined above in (|4.19|) . 

Acting on (^2l]) - (|l2^) with 8^ for £ = 0,1,2 and using: 



d* i L k (<f> i ) = {ML*) (L k (^)) + [k+^)L 



-k+e-i 



(*0 



(4.24) 
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valid for £ = 0, 1, 2 and k > 0, we get: 

'ml' 



31 X, 



(i) 



Mi 



v(o) 



+ kX. 



(o) 



+ (fc - £ + 1) X 



9SX; 



(2) 



(ml') , X 



(2) 
k 



(1) 

(2) 



+ {k-2{£-l))X^> 



If -(£-!)) Xi 



.(o) 

L fe+£-i 



(4.25) 
(4.26) 

(4.27) 



Here we recognize a structure of Wi +00 algebra under substitution I — > £ — 1. Now we are ready 
to restrict attention to the part of the algebra involving A^j^ operator from ( 4.22j ). We note that 
( 4.18] ) and identity (A.4) from Appendix [A] enable us to obtain: 



vi 1 ) T 
A k-i ' L 



k m 



(4.28) 



i=l 



Our main result is contained in the following: 



Proposition 4.1 The correct additional-symmetry flows for the cKP rm hierarchies (2.31), span- 
ning the Virasoro algebra, are given by: 



dZL 



i.e., with the following isomorphism 



ML k ) +X£} 1 ,L 



(4.29) 



+X£} V where X^\ are defined in (gj||). 



Accordingly, on dressing operators and BA functions the flows (4-2&) read: 



dtW 



(ML k ) _+X^ W ; flj VM) = (" (ML k ) _+X^ (^(i,A)) (4.30) 



Indeed, using ( 4.28j ) we first find that (d k f L)_ can be cast in the form of ( [4.12 ) with: 

= (ML fe ) + ($ i ) + ^ fc - 1 (*i) + 4 1 i(*i) 



K ),(*<) + « 



*\k— 1 , 



(1) 
fc-1 



(4.31) 



Taking into account that X\\ = for i = 0, 1, 2 we see that eq.( 4.31 ) reproduces ( |4.13|) - (f4.15| ) 
(with ambiguity on the right hand side of (4.14) removed by fixing a = /3 = 1/2). Hence d\ = di\ 
for £ = 0,1,2. 

Secondly, we note that the modified additional symmetry flows defined by (4.29) commute with 
the isospectral flows (2.3) and, due to ( 4.31 ), they preserve the form of tfie cKP rim Lax operator 
(2.31). The remaining question is whether they form a closed algebra. Indeed, using identity (4.26) 
we arrive at the fundamental commutation relations for £ = 0, 1, 2 and any k > : 



[d* £ ,d* k ]L = (k-£)d* k+e _ 1 L 



(4.32) 
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This discussion shows that [£% , Ck] = (i — k)£i+k for i = —1,0,1 (si (2) generators) and 
arbitrary k, where Ck-i ~ —d^ . Since according to proposition 11 the generator £2 is associated 



with X2 — (ML 3 )_ , all higher Virasoro operators can be obtained recursively from: 

Cn+i = - t ~* [A., A] , n>2 (4.33) 
(n-1) 

Then eq.( (4.32|) implies that £ n with n > 3 may differ from the generators given by the flows 
d^n+i ~ ~~ (ML n+1 ) _ + Xn^ defined in ( 4.29 ) at most by flows commuting with the sl(2) additional 



symmetry generators, i.e., by ordinary isospectral flows. Therefore, we can now easily show by 
induction that £k , k > —1, obtained in the above way form a closed Virasoro algebra up to 
irrelevant terms containing ordinary isospectral flows. 

For illustration consider, e.g., [£2 , £3 ] = Z. Commuting with Z we find [£-\ , Z] = 6£^, 
which fixes Z to be —£5 up to isospectral flows commuting with the sl(2) subalgebra. It is easy to 
see how to extend this argument to cover the whole Virasoro algebra. 

4.3 Additional Symmetries versus DB Transformations for cKP r m Hierarchies. 
String Condition 

Let $ be an eigenfunction of L defining a Darboux-Backlund transformation, i.e. : 

_ $ = L^($) , L = (^D®- 1 ^) L (^D- 1 ®- 1 ^ , W = (^D®- 1 ^) WD' 1 (4.34) 



Then the DB-transformed M operator (cf. ( |4.4| )) acquires the form: 

M = (^D^ 1 ) M (fcir 1 ^- 1 ) = J2 l -^t r+l l\_ + M_ (4.35) 



l>0 



M_ = WX^W- 1 -t r -J2 —tr+i-i-W . W- 1 (4.36) 

l>i r atl 



where Xr r \ = DX^D 1 with X^ as in (4.3). Clearly Xr r \, like Xr r \, is also admissible as 



canonically conjugated to D r . The DB-transformed M-operator acts on the DB-transformed BA 
function ( 3.1 0|) as: 



Af^(A) = ( — + 5r (A) ) tP(X) , 5 r (A) = a r (X) + -A" r (4.37) 
\ r OA I r 



Now, let again L belong to a cKP rjm hierarchy (2.31) and let us consider auto-Backlund trans- 



formations djUHHH) preserving the constrained cKP rm form of L : 



L = T a LT- 1 = L + + Y / ^iD- 1 ^ l , (4.38) 
i=i 

$ a = T a L($ a ) z ^a = <S>a 1 (4-39) 

^i = T a (^) , * f = 37 1 ** i = -*- 1 5- 1 (* i * a ) , i + a (4.40) 

where the DB-generating <3? = <£ a coincides with one of the eigenfunctions of the initial L fl2.31|) . 
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With the help of identities ( |A.5|) -( [A"lj| ) from the Appendix |A] we find the following explicit form 
of the DB transformation of the operators -Xl_i ( 4,22 ) : 

TaX&Ta 1 = ^-(£W)^V{T 8 (x« 1 + ^*- 1 )(* a )}D- 1 *J 1 (4.41) 



fc-i 



fc-2 



J=0 



(4.42) 



Here L, T a are as in ( |4.38| ) and the DB-transformed Xu~\ have the same form as X£± x in (gjg) 
with all (adjoint) eigenfunctions substituted with their DB-transformed counterparts as in ( [4.38| )- 

( 4.40| ). Also notice that in the particular case of cKP r i hierarchies | LW) ( p~42|) coincides with 

the (pseudo-differential part of the power of the) full cKP r i Lax operator (cf. eq.( 2.31| ) for m = 1 
and jAlD ). 

Taking into account ( f4.38| )-( |4.40| ) and ( 4.41| )-( [4.42| ) we arrive at the following important 

Proposition 4.2 The additional- symmetry flows ( \j-2!\ ) for cKP r j hierarchies (eg. (ft. 31 ) withm = 
1) commute with the Darboux-Bdcklund transformations ( £.38 ) preserving the form of cKP r i, up 
to shifting of fy.2i> ) by ordinary isospectral flows. Explicitly we have: 

dl 



-(I) 



diL 



(ML k ) 



+ 



dt k - 



(4.43) 



Proposition L2 shows that the additional-symmetry flows ( [4.29| ) are well-defined for all cKP r] i 
Lax operators belonging to a given DB orbit of successive DB transformations. Notice that it is 
precisely the class of (truncated) cKP rj i hierarchies which is relevant for the description of discrete 
(multi-)matrix models (refs.(2^, |24|, fl8| and sections below). 

Motivated by applications to (multi-)matrix models (see section ||), one can require invariance 
of cKP rjm hierarchies under some of the additional-symmetry flows, e.g., under the lowest one 
<9q = <9o,i, known as "string-equation" constraint (string condition) in the context of the (multi- 
)matrix models: 



d* n L = 



Mi , L 



-1 



<9n$ = 



M + $ = 



(4.44) 



Eqs.( 4.44 ), using second eq.( |4.2|) ,(4.4) and first eq. (|4.31| ) for k = 0, lead to the following constraints 
on L ( |2.31[ ), the BA function ij}(t,X) and its DB-generating eigenfunction &(t), respectively: 



El + r d 

1>1 r dtl 



^l + r d 




-1 

A 1-r d 
r d\ 



(4.45) 
(4.46) 
(4.47) 



Recall now formula ( 3.54 ) for the r- function of the cKP rim hierarchy ( 2.31 ). Noticing that the 
eigenfunctions <3?( fc ) of the DB-transformed Lax operators satisfy the same constraint eq.( 4.47j ) 
irrespective of the DB-step k, we arrive at the following result ("string-equation" constraint on the 
r-functions) : 
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Proposition 4.3 The Wronskian r -functions l \3. 54 ) o/cKP rm hierarchies (2.31), invariant under 
the lowest additional symmetry flow ( 4-44 )> satisfy the constraint equation: 



^l + r d \ r( n ) 



n-l 



(0)1 



(4.48) 



In particular, eq. ([4.47 ) applied to (|3.32|) , (|3.33| ) leads to the following restriction on the Darboux- 
Backlund generating eigenfunction $j° (t) of the "initial" which uniquely fixes the form of the 
members of the pertinent cKP rm DB-orbit: 



Corollary 4.1 The constraint ( 4-4 r > ) imposes the following relations on the eigenf unctions defining 
the cKP rm DB-orbit: 



4>i^ (A) = const 



4 k) (t) = f dXX k ^^(t,X) = f dAA Kl+O r(fc) ^ - [A ^ ^(t.A) 



r( fe )(t) 



(4.49) 
(4.50) 



$f)(t) = / dXX k ^ k) (t,X) , i = 2,... 



m 



5 One-Matrix model, DB- Wronskian Technique and Virasoro con- 
straints 



The standard (Hermitian) one-matrix model has partition function |^9[ : 

Z N [t] = J dMexpljrt r TrM r \ (5.1) 



,r=l 



where M is a Hermitian N x N matrix. One uses the method of generalized orthogonal polynomials 
J60| ] to evaluate partition function fl5,l| ). After angular integration in ([5[l]) we obtain (up to an 
overall constant) : 



- / n^exp n - A i) 2 ( 5 - 2 ) 

i=l l \i=l k=l I ) i>j=l 



z N [t] m 

Comparing ( |5.2| ) with (3.89), we notice that (5.1) precisely coincides with the r-function of 
the generalized BH hierarchy, i.e., the Darboux-Backlund orbit of cKPn hierarchy starting from 
a "free" Lax = D, with the special "initial" condition fl4.49| ) : 

0<°)(A) = const (5.3) 

which, as discussed in section |] reflects the imposition of the additional requirement of invariance 
of the cKPi i solutions w.r.t. the lowest additional symmetry flow. 
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Reversing the argument leading from to ( |5.2[ ), we find that the r- function of a generic 

cKPn Darboux-Backlund orbit (eq.( |3,47l) for r = 1), i.e., a Darboux-Backlund orbit starting from 
a non-free initial Lax = D + <J>(°).D _1 \iir( ) anc l with <E>(°) having an arbitrary "density" ^>( )(A) 
not subject to the constraint (|5.3|), can be identified with the partition function of the following 
generalized one-matrix model: 



Z N [t] = J dMexp|Trln0 (o) (M) + £t r .TrM r + lnr (0) (i-Tr[M" 1 ])| 

r (0) (t -Tr[M^]) = r (0) (ti - Tr M~ 1 ,t 2 - -TrM~ 2 ,£ 3 — - Tr M~ 3 , . . .) 

2 3 



(5.4) 



^lnr(°) = cD( )M/( ) 

Below we will exhibit in more detail the relevant integrability structure leading to the Darboux- 



Backlund Wronskian solution ( |3.8E ) for ( |5,1| ) and its generalizations (| 

One deals with the integral (|5.2 ) (or the corresponding generalization from (|5.4j )) using orthog- 
onal polynomials: 

P n {\) = X n + o(x n ~ l ) , n = 0,l,... (5.5) 
which enter the orthogonality relation: 

J d\P n (X)P m (X) exp (j£ tk^J = h m (t)5 nm (5.6) 



From ( |5.5[) and (5J3) follow the recursion relations for the orthogonal polynomials: 

oo 

XP n (X)=J2QmPi(X) (5.7) 

1=0 



The matrix elements of Q can be found by comparing (|5.7| ) with the definition of orthogonal 
polynomials, yielding: 

Qn,n+1 = 1 ; Qn,l = / > n + 2 ; Q n>n -l = 71 > 1 (5.8) 

The orthogonal polynomials approach leads to expression for the partition function in terms of h n : 

Zn = const hohi ■ ■ ■ /ijv-2^v-i (5-9) 



From the orthogonal relation ( |5.6| ) and definitions (5.5) we also obtain: 

<91n h r , 



dt r 

as well as: 



(Qlnn (5-10) 



f)P n_1 

^r = -zZQ r m p iW (5.11) 

olr 1=0 



Defining now a wave function: 



y n (t,X) = P n (X)e v ^ , V(t,X)^Y, xht k ( 5 - 12 ) 

fc>i 
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we obtain from (5.7) and (5.11) the following eigenvalue problem: 

A<£ = Q^> ; 



— = Q\ + )* 



Of, 



(5.13) 



where ^ is a semi-infinite column (\^0) ■ • • ; ^ni ^n+l> ■ ■ -) T an d Q is a semi-infinite matrix i.e., with 
indices running from to oo. Furthermore we adhere to the following notation: the subscripts — /+ 
denote lower/upper triangular parts of the matrix, whereas (+)/(—) denote upper/lower triangular 
plus diagonal parts. The compatibility of the eigenvalue problem ( |5.13D gives rise to a discrete 
linear system: 

d 



at 



-Q — Q 



(5.14) 



Choosing parametrization Q nn = ao(n) and Q n ,n-i = o>i{n) = for the matrix Q, we can 

rewrite ( 5.13j ) as: 



AW n = ^ n+ i + a (n)* n + ai(n)^ n _i 
d^ n = + a (n)^ n (for r = 1) 



(5.15) 
(5.16) 



Consistency of ( 5.16 ) and ( [5,15 ) yields the following evolution equations having form of the Tods 
lattice equations of motion: 



dao(n) = a\{n + 1) — ai(n) 
da\{n) = a\{n) (ao(n) — ao(n — 1)) 

Using an inverted version of ( |5l6l) : 



(5.17) 
(5.18) 



(5.19) 



d - a (n) 

one can re-express the linear problem ( |5.15[ H( |5,16| ) entirely in terms of \l/ n at a fixed lattice site n 



_d_ 

8t r 



Af n = L n ^ r 
*n = (L n )[ 



n > 



where now L n is a continuum pseudo-differential Lax operator associated with the site n : 

1 



L n = d + ai(n)- 



(5.20) 
(5.21) 

(5.22) 



d — ao(n — 1) 

Note that from ( 5.16| )-( |5.15| ) we also have X^ n +i = ^n+i^n+i) J '- e -i the lattice shift n — > n + 1 



does not change the eigenvalue A. In fact, the lattice shift n — > n + 1 can be given a meaning of 
a Darboux-Backlund transformation within the generalized Burgers-Hopf hierarchy introduced in 
section [3|. To see this, we rewrite ( 5.19| ) as follows: 

* n+1 = e f ao ^d e-/ ao(n) ^ n = ^(n)^- 1 ^)^ = T{n)^ n (5.23) 

where 3>(n) = e-f a °( n \ and T(n) = < I>(n)5 < I ) ~ 1 (n) plays a role of the DB transformation operator 
generating the lattice translation n — * n + 1. Combining ( |5.23j ) with ( |5.15 )-(5.20) we find: 



L n+ i^ n+ i = \(d - a (n))^ n = (d - a (n)) L n ^ r 
= (5-a (n))L ri (9-a (n))" 1 * ri+ i 



(5.24) 
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Thus, the Lax operators at different sites are related by a DB transformation: 

L n+1 = (8 - a (n)) L n (d - ao(n))- 1 = T^L^in) 
with the Lax operator ( |5.22 ) taking the form: 



L,, 



$(n) 



L(n) = d+ $(n)0 -1 tt(n) 

e / a o(n) =ai ( n)e /a (n-l) > ^ _ - J a (n-l) = ($(n _ 1 



(5.25) 

(5.26) 
(5.27) 



This is precisely the form of the Lax operators belonging to the cKP^i DB orbit starting from a 
"free" Lax Lq = d ( 3.87 )-( 3.88 ) or, according to subsection 3.3, L n ( |5.26|) describes the generalized 
Burgers-Hopf system. 

We see that we proved: 



Proposition 5.1 The one-matrix model problem (5A) is equivalent to solving the generalized 
Burgers-Hopf system. 

In more detail, realizing from ( [5 . 1 Q| ) that ao(n) = dlnh n and accounting for ( [5,12 ), we have the 
following identification between discrete one-matrix model objects h n , P n (A), Zn , on one hand, 
and eigenfunction $(n) = <3>( n )(t), BA function i/j^(t,X) and r-function of continuum n-truncated 
KP hierarchy (Prop,|3~7^), on the other hand: 



h r . 



<£(n) 



tf n (t,X) = P n (X)e v ^ = X n ^ n \t,X) = (^wf ] {t)X n -iy 

3=0 

Z N = hoh 1 ---h N - 2 h N -i = *(0)---$(JV-l) =W N [h ,dh ,...,d N - 1 h } (5.30) 

In particular, from ( |5.29D and ( |5.30 ) one obtains an explicit expression for the coefficients of the 
orthogonal polynomials (cf. refs.(6l], |62], p4|1): 



(5.28) 
(5.29) 



PnW 



X 



V x r 

p W n [h ,...,d^h ] 

n W n [h (t - [A- 1 ]), dh (t - [A- 1 ]), . . .,d n - l h (t - [A- 1 ])] 



W n [ho(t),dho(t),...,^h (t)} 



Further, using an identity: 



= | dX-^ jexp (Jr P n (A)P m (A) j 



we derive: 



~dX 



P* ; P=J2kt k Q k + - 1 

k=2 

Consistency with ( |5.13 ) requires the so-called discrete string equation: 



[Q,P] 



oo 
k=2 



Of 



k-l 



(5.31) 



(5.32) 



(5.33) 



(5.34) 
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which using ( |5.10| ) can be rewritten as: 



2^kt k - hxh n =0 n>0 (5.35) 



dh dtk-i 



According to ( |4.47| ) and the identification ( p\28|) , eq.( |5.35| ) is the discrete one-matrix model coun- 
terpart of the "initial-condition" constraint ( |5.3j ) on the associated Burgers-Hopf hierarchy due to 
invariance under the lowest additional symmetry flow. 



In refs.[19|, using Ward identities' techniques for the matrix integral ( |5.1|) , it was shown that 
the partition function Zn in fact satisfies an infinite number of constraints spanning the Borel 
subalgebra of the Virasoro algebra: 

C { S N) Z N = , s > -1 (5.36) 

4~> = g*£ + W £ + |A_»_ , s >, (5 ,T) 

Here we provide a simple alternative proof that the Wronskian solution ( 5.3C| ) indeed satisfies all 
the Virasoro constraints: 

Proposition 5.2 The Wronskian T-junction of the generalized Burgers-Hopf hierarchy §3.890 sub- 
ject to the "initial" condition A5J\) satisfy the Virasoro constraints: 



£i N) W N [<j), d<f>, ... , d N ~ l (l)\ =0 , s > -1 (5.39) 



Proof. The case s = — 1 in (5.39) is a particular case of Prop. 4.3. For s > the proof will go 
by induction w.r.t. N and we shall employ the r-function recurrence relation fl3,43|) which in the 



special case of (3.8S) takes the form: 



W N +i[4>,d4>, ... ,d 



N 



d\\ 2N :e- § ^ :W N [cf>,< 



■ 



N-l 



(5.40) 



Using the explicit form of the vertex operator with 0(A) as in ( 3.44 ), it is easy to deduce the 
following commutation relation: 



£i Y - J 1 : <- f),X> :-: r"^' : 1 I ( A* +1 ^ + (2N + s + 1)A«) : e"* A > : 

,d_ 
dt k 



+ :e- § ^ : A* £ A" 



k=l 



Applying ( |5.41 ) on both sides of ( 5.40 ) and integrating by parts we get: 



4 JV+1 %^ +1 [<A,^,...,a 



N 



d\\ 2N :e- § W :CWw N [<f>, 



N-l 







+ /dAA^: e -^):(]TA- fc ^-)T^v[ 



N-l. 



k=l 



(5.41) 



(5.42) 
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Thus, to complete the proof we have to verify that the last term on the right hand side of (5.42) 
does vanish identically. Indeed, this term can be rewritten using second equality ( 3.8SD as: 



N 



N-l 



N 



N 



n ^ e e (^r k n (a, - a,) n w n 

j=0 k=l 1=0 i>j j=0 k=0 



fc )e«*.^)} (5.43) 



valid even without the restriction ( |5.3| ). The integral (5.43) vanishes due to the fact that the 
expression in its integrand: 



s N-l N 

E E ( A o fc (Aivr k n ^ - A i) n M 

k=l 1=0 i>j j=0 



(5.44) 



is anti- symmetric w.r.t. at least one pair of A's. □ 



6 Two-Matrix Model and cKP r=P2 i m= i Hierarchy 

6.1 Two-Matrix Model and the Corresponding Linear Toda-like System 



We consider the two-matrix model with a partition function [50 : 

Pi P2 



Z N [t, t, g}= J dMxdM 2 exp J ]T *r Tr M i + E ^ Tr M 2 + 9 Tr M i M 2 | (6- 1 ) 



^r=l s=l 



where are Hermitian N x N matrices, and the orders of the matrix "potentials" pi t 2 may be 
finite or infinite. After angular integration in (|6.lD one obtains |30| (up to an overall constant) : 



1 r { 1 

ZN[t,t,g] = — n dA idA i expiE( y (Ai) + ^(Ai)+5A i A,) > A(Aj)A(Aj) (6.2) 

i=i yi=i ) 

where V(A) = EfcLi*fcA fe , V(X) = Esli^A 8 . A(Aj), A(Aj) are standard Van der Monde determi- 
nants. 

Similarly to the one- matrix case, we can write down an explicit Wronskian expression for ( |6.1[ ). 
To this end we notice that the integral over one set of eigenvalues, e.g. Aj, can be represented as: 



N N 



I IT II (A* - h) II {p( x h, *k)e^>~ Xk) } = W N [$(t; Ai), . . . , Ajv)] (6.3) 
i=i i>j k=i 

$(t;A) = f dXp(X,X)e^'~ x) , p(A,A) = e sAX ; x = f i (6.4) 



Here and below £(t, A) = E«=i ^A s and A) = ErLi *rA r . Inserting (|o|) into ( p!2[) and using 
the simple identity ( 3.49[ ) we obtain: 



Z " = / ft dXi IT (Ai - Ai) W N [e^^(t; AO, ... , e«'' A *)*(i; Ajv)] 

i=l i>j 
r N N 

= Jl[dX i ]l(X 3 y- 1 W N [e^^t;X 1 ),...,e^^(t-,X N , 

8=1 jf = l 

a , d N ~ 1 
oh dh 



W N [ho, —ho, - jV _ 1 ho] (6.5) 
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where: 

ho = h (t, t) = I dX dX p{\, A)e ?( *' A)+f ( '- A) (6.6) 



Since in the derivation of (|6.5| ) the "density" p(X, A) was in fact arbitrary (not necessarily fixed 
by the expression in (|6.4D), eq.(|6.5|) yields simultaneously the Wronskian solutions of all higher 



multi-matrix models |19| provided we choose p(X,X) in the form of eq.(7.9) from next section. 



In ref.[31] it was shown that, by using the method of generalized orthogonal polynomials pO 



the partition function (|6.l| ) and its derivatives w.r.t. the parameters (t r ,t s ,g) can be explicitly 
expressed in terms of solutions to generalized Toda-like lattice systems associated with ( p\l| ) . Unlike 
the case of the full generalized Toda lattice hierarchy ||], however, the associated Toda matrices 
are (i) semiinfinite, and (ii) contain in general finite number of non-zero diagonals. In the present 
section we present in detail the identification of the correct continuum integrable system underlying 
the explicit Wronskian solution ([3J^) of discrete two-matrix models in the case when at least one 
of the matrix potentials is of finite order, e.g., P2 = finite. 

As in one-matrix model one can deal with the integral fl6.2|) using generalized orthogonal poly- 
nomials [60] : 



P n (X) = X n + [ A ; P m {X 2 )=X m + A" 1 " 1 , n,m = 0,l,... (6.7) 



which are orthogonal to each other with respect to the weight exp (^U(A) + V(X) + gXXj : 

h n 5nm = J fdXdX P n {X) exp I ^ A% + £ X s t s + 5 AA j P m {X) (6.8) 

Here h n are the normalization factors. In this setting of generalized orthogonal polynomials the 
multiplication by A, A can be represented by the Jacobi matrices Q and Q: 

n+l n+1 

xp n {x) = Y,QnM x ) ; xP n (x) = Y J PiWQiAhn/h l ) (6.9) 
1=0 1=0 

Qnh Qui are matrix elements of Q and Q, which define in what follows the underlying linear system. 
Due to definitions Q6.7| ) they satisfy the relations: 

Qn,n+i = 1 ; Qn,i = / > n + 2 and Q n +i,n = h n+ i/h n ; Qi >n = l>n + 2 (6.10) 

The orthogonal polynomial approach leads to expression for the partition function in terms of 
h n which is similar to the one-matrix case (|57 



iV-l 

Zjy = const | [ h n (6-H) 

n=0 



From the orthogonal relation ( |6.8j ) and the recursive relations ( |6.9| ) we obtain: 
and, accordingly: 

8 Af-l o N-l 

w \nZ N =J2 Qn,n , gf^Z N =J2 Qn,n (^) 
r n=0 s n=0 
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In terms of matrices Q and Q the corresponding linear problem for the lattice wave function 



ip n = P n (X) and the Lax (or "zero-curvature") representation of the latter read [31] : 

8 l 



(6.14) 

ip m (6.15) 

n 

(6.16) 
(6.17) 

where the Toda matrices are subject to the following additional constraints (known as "coupling 
conditions" [p]] or discrete "string equation" ) : 



n, m 1pm = A-^n 




1 






-(QL 


)nm ^ m 


• at*" 


dt r ® 


Q\+) , 


Q_ 


? 


at. 


dt r ® 


Q\+), 


Q 


? 


at. 



1 

Qh 

Q(+) 



-g 



[Q,Q 



V 1(s + 1) t o 

8=1 
Pi — 1 



1~ -n 
-till 

5 



^ - — - — -*r+lQ(+) — — ^ll 



<7 



r=l ^ 

The parametrization of the matrices Q and Q is as follows: 

<3n,n = ao(n) , Qn,n+1 = 1 , Qn,n-k = Ofe(n) fe = 1, — , jP2 — 1 
Qn,m = for m — n > 2 , n — m > p2 
Qn,n = b (n) , Q n ,n-1 = Rn , Qn,n+k = M^O^n-fl ' ' ' "^n+fc k = 1, . . . ,p 1 - 1 

Qnm = for n — m>2 , m — n>p\ 



(6.18) 
(6.19) 

(6.20) 



(6.21) 



(6.22) 



The Q, Q-matrix elements can in general be expressed in terms of the normalization factors h n 



from ( p.8| ). For some of the elements the relation to h n and Z n takes a very simple form [31] : 



a (n) 



d 

lnh n , b (n) = -~- lnh n , R r , 
ot\ 

d — i q N—i 

— lnZ N = V ao(n) , -^lnZ/vr = V &o(n) 
ail „ oh ^~ ' 

n=0 



JV-1 



/?.„ 



*n-l 



(6.23) 
(6.24) 



n=0 



We now exhibit the Toda structure behind the linear problem ( |6.14 )-( |6.17| ) focusing on the lowest 
flows with r,s = 1. We start with the second lattice equation of motion (|6.17| ) for s = 1. Using 
parametrization ( |6.22 ) we obtain: 



- r R n = R n (b (n)-b (n-l)) , b Q {n) = h{n) - h{n - 1) (6.25) 
ot\ at i 







h(n) 



dt\ \R n +i • ■ ■ R n +h 



b k+ i(n) - b k+1 (n - 1) 

Rn+l ■ ' ' Rn+k 



k > 1 



(6.26) 



Similarly, the first lattice equation of motion ( 6.17| ) for r = 1 gives : 



d d 
— b (n) = R n+ i - R n , —bk(n) = R n+ ib k -i(n + 1) - R n+k b k -i(n) (6.27) 
ot\ ot\ 
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for k > 1. In complete analogy, the lattice equations of motion ( 6.16) ) for r = l,s = 1 read 
explicitly: 



d d 

— r a Q (n) = R n +i-Rn , T^a k (n) = R n _ k+1 a k -i(n) - R n a k ^i(n - 1) 

Oti Ot\ 



(with k > 1) and: 
9 



ao(n) = ai(n + 1) — ai(n) 



a fc+ i(n + 1) - a fc+ i(n) 



<9ii 5ti \R n - ■ ■ R n - k+ i J R n ■ ■ ■ R n - k+ \ 

with fc > 1. Following ( ggg ), ( pT27p , flOSp we obtain the "duality" relations: 



5 

— 6i(n) = 
oil ail 



d d d d 

^r a o(n) = ~r&o(") , ^r«i(«) = irr^n 
ail oti oti oil 



From the above one gets the two-dimensional Toda lattice equation: 
d d d 

— InRn = a (n) - a (n- 1) -> — - r lnR n = R n+ i-2R n + R n _i 
oti Ot\ oti 



Also, using (|6T29|) , (|Q^) and ( fT27D , we get from fu|) : 



9 2 



In Zn 



dtidt 



N 



ai(N) , 
Rn 



2 



lnZ N = h(N-l) 



dh 



h N - 



7 2 

N 



(6.28) 



(6.29) 



(6.30) 



(6.31) 



(6.32) 



We note the two-dimensional Toda- like lattice structure exhibited by equation ( |6.32| ). 

For later convenience, let us also display the lattice equations of motion for powers of Q and Q 



d r 
— Q r = [Q (+) , Q r 

d 



Q s — Q(+),Q S 



dh 



as well as for their matrix elements 



4-Q Q'.Q- 

Oti 

d 



dit 



Q s = Q S ,Q- 



—Q s 



_d_ 

dh 



Q 



n,n—k 



dh Wnn 



Q S n. 







dh \ R n • • • R 



n—k 



Rn+lQ n ,n+l ~ ^nQn-l.n 

Rn— k+lQn,n— fe+1 ~~ R"nQ n —l,n—k 

Qn+l,n ~ Qn,n-l 
Qn+l,n-k ~ Qn,n-k-l 



k = 1, . . . , S 



TO — fe+1 



R n - ■ ■ R 

n— fe+1 



k = 1, . . . , s 



(6.33) 
(6.34) 



(6.35) 



(6.36) 



Let us note that, due to the upper Jacobian form of Q having only one non-zero lower diagonal, 
the matrices Q s have s non-zero lower diagonals, i.e. 



Qn,n- 



for k > s + 1 
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s-1 

Qn : n-s = Rn ■ ■ ■ R n -{s-l) , Qn,n~(s-1) = Rn ' ' ' R n -(s-2) ^ b (n - j) (6.37) 

j=0 

Qn,n-(s-2) = Rn---Rn~(s-3)\J2 b ^ n -3)+ Yl b {n - i)b {n - j) 

\j=0 0<i<j<s-2 

For completeness, let us also add the explicit expressions for the higher flow equations for bo(n 

d 



bi(n), Rn+i resulting from ( 6.16 ) and ( 6.17 ) : 

d /~«\ 



T^bo(n) 



_d_ 

dts 




Rn+l 



W (Q S ) , - r b l {n) = - r [R n+l (Q< 
Oti v 'nn dt s Ot\ \ v 

l(Q s ) 

t 1 Jn+l,n 



n,n+l 



_0 

at 
o 



d_ 

dt r 



R 



n+l 



d_ 



+l,n 



(6.38) 



(6.39) 



(6.40) 



which easily follow by comparing equations ( 6.16 )-( 6.17 ) with equations (|6.33|) - (|6.34j) , respectively, 
for the main diagonal and the first upper /lower diagonal matrix elements. 

6.2 Equivalent Hierarchies Method and the String Equation 

The notion of equivalent hierarchies was introduced to answer a question of equivalence between 
different formulations of integrable models. The key role in this context is played by the particular 
form of coordinate transformations on the space of evolution parameters: 



t n —> t n — t n + O (£n+l) tn+2, ■ ■ •) 



(6.41) 



This coordinate transformation appeared in [45] as a tool used to establish an equivalence between 
Zakharov-Shabat hierarchy associated to the zero-curvature equation : 



dB n 3B r , 



dt r , 



Otr 



+ [B n ,B m ]=0, 



n,m = 1,2,... 



(6.42) 



and the standard KP hierarchy formulated in terms of pseudo-differential operators. Because 
the zero-curvature equation ( 6. 42] ) is coordinate free, the differential operators B m defined by the 
transformation ( |6.41D : 

Otr 



Bn(t)=Y,^B m (t) 



(6.43) 



will continue to satisfy equation ( |6.42| ). As shown in [45], there exists a unique transformation of 
the form (6.41) for which the differential operator B m takes a standard form of (L)™ with the Lax 
operator L = D + u %B~ l being a solution of the Lax equations of the KP hierarchy with respect 
to times i n , n > 1 : 

81 

= [B m ,L] , B m = (L)™ (6.44) 

Otm, 



The general conclusion is that the invariance of the hierarchy under ( 6.41 ) can be used to cast 



the differential operator hierarchy defined by (6.42) into the particularly convenient standard form 



43 



of the KP hierarchy. Similar considerations have also been made in 45] about the Toda lattice 
hierarchy. 

The method of equivalent hierarchies appeared also in the setting of the Generalized Kontsevich 
Model pi. 



In ref.[29] a similar in spirit notion was proposed to deal with the string equation constraint in 
the setting of the two- matrix model as represented by the Toda system ( |6 . 1 4| ) - ( |B . 1 T| ) augmented 
by ( 6.18[) -( |6.20 ). In this reference an equivalent Toda hierarchy was found for which the string 
condition could be partially implemented by a simple identification of two flows. We now expand 
further on these ideas. 

The first step is to introduce the fractional powers of Q. 

Definition 6.1 Matrix Q is defined as the fractional power of the Q-matrix according to: 

Q = Q^r (6.45) 



Let us note that it is at this point where the finiteness of the order p2 of the second matrix potential 
in ( |6.1[ ) becomes crucial. 

The matrix Q is, unlike Q, an upper Jacobi matrix with only one non-zero lower diagonal, 
therefore, Q has the same matrix form as Q and it is natural to parametrize it in the analogous 
way (cf. eq.flo^)) : 

Qnn = bo(n) , Q n ,n-1 = Rn , Qn,n+k = bk(n)R n +± ■ ■ ■ R n ^ k k = 1, . . . , p\ — 1 

Q nm = for n - m > 2 (6.46) 

Correspondingly, the matrix elements Q s n n _ k have the same expressions as given by ( |6.37[ ) with 
b (n), bi(n),R n replaced by b (n), bi (n), R n . 

Furthermore, the lattice equations of motion for powers of Q: Q s = Q^- 1 ; which follow from 
( 6,16] ), take the form: 

d ^ s 

~g£~Qnn = Rn+lQ n ,n+l ~ RnQn-l,n 

d ~ s 

T£~Qn,n-k = Rn~k+lQn,n-k+l ~ RnQn-l,n-k k = 1, . . . 8 (6.47) 



{ft Qnn Qn+l,n Qn,n—1 



d_ 

9 I Qn t n-k \ Q n +i t n-k ~ Q n ,n-k-l 

dt\ \ R n ■ ■ ■ R n -k+l I Rn' ■ ■ Rn-k+l 



k = l,...,s (6.48) 



The system ( p.47[ )-( |6.48 ) is exactly of the same form as the corresponding equations for Q ( |6.35| )- 

Starting with eqs.( 6.47| )- (|6.48; ) for k = s and using properties ( 6.37| ) (and their analogs for Q), 
we arrive by induction w.r.t. k = s, s — 1, . . . , at the following result: 

Proposition 6.1 The matrix elements of Q are completely expressed in terms of the matrix ele- 
ments of Q through the relations: 

= E7s CT Qf_) , a = 0,l,...,p2 (6-49) 

cr=0 
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where the coefficients 7 Sj o are ti-independent, whereas the coefficients 7 scr with a > 1 are indepen- 
dent ofti,t\. All 7 S(J are simply expressed through the subset thereof: 

s + l~ , . 

a s = 7p2-i,a = — t s +i (6.50) 



The equality ( 6.50| ) follows from the coincidence of eq.(6.49) for s = P2 — 1 with the "string equation" 

(Hi). 

Let us note that, although the system of equations ( 6.4S| ) for Q is overdetermined, it is com- 
patible by construction. A calculation yields: 



7ss = (7n) s , ls,s-i = s (7n) s Sio , 7n 



P2r \P2-i 

„ P2 



(6.51) 



ls,s-2 = (711 
731 



■ s-2 



s ( s - !) ,2 



37ii 



7io 



2 



(710) + s 
1 



1 37 710 



11 



, 710 



■ £p2-l ( ' 



P2~ 2 

P2r \ p 2-i 



"P2 



5(P2 - 1) 



^2-2 



P2 ~ 1 ^a-l 
2^2 L 2 



f — Elf Va -1 

Relation ( |6.49| ) implies the following expressions for the Q = Q^- 1 -matrix elements: 



(6.52) 



R n = JnRn , b (n) = 711 6 (n) + 710 , h(n) = 711^1 ( n ) + 



731 

3711 



7io 



(6.53) 



etc.. 



The preceding discussion can be rewritten in the dual form with the roles of Q and Q inter- 
changed. In particular, the matrix elements of Q are completely expressed in terms of the matrix 
elements of Q through the relations: 



r = 0, l,...,pi 



(6.54) 



p=0 



Here, the coefficients rj r fi are ti-independent, whereas the coefficients r\ rp with p > 1 are independent 
ofti,ti. All expressed through the subset thereof (cf. eq.( |6.20D ) : 



0r — Vpi — l, 



r + 1 



W+l 



(6.55) 



in the same form as ( 6.52] ) with j S(7 , a s , s,p2 replaced by rj rp , f3 r ,r,pi, respectively. 

Combining relations ( |6.49| ) and ( 6.54| ) lead to the following explicit expression of Q in terms of 
Q or vice versa: 



Q 
Q 



a 



P2-l ( 



1 Pl-1 
r=0 

1 P2-1 _ 

^-i^): 1 " 1 + e ^qu 

s=0 



(6.56) 
(6.57) 



with a s ,l3 r as in fl6.50| ),( 6.55 ). 

The next step, natural from the viewpoint of equivalent hierarchies (see equation ( |6.41| )), in- 
volves introducing a new subset of evolution parameters instead of {t s }- 
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Definition 6.2 The new subset of evolution parameters j^sj * s defined through the relations: 

(6.58) 



9 -ST 9 



which imply: 



(6.59) 



t s — t s (isi is+i 5 • • • j 
i.e., t s does not depend on t\, . . . , t s _i. 

In particular, from ( |6,51 )-( |6.52| ) we have a number of technical relations (to be used later on) : 

d d -\ ( Vix p 2 

= 7u^r > 7n = — t P2 \ =— — t P2 

dh dh \ 9 J 9{P2 ~ 1) 

(6.61) 



(6.60) 



7io 



>p 2 



(p 2 ~ l)t p2 -l 731 2 _ (P2 - 2)tp 2 -2 

7io — ~ 



P2tp 2 
9 ( P2 



-u 



3711 

-(P2-1) 



P2t 



in 



P2 \g(p2 - i)" p - 

-(P2-2) 



P2-1 



P2 



5(^2 - 1 



P2 



"(P2-1) 



^P2 — 1 



P2 



~t)i 



tp2-2 



P2 , pa-l ^ 



.5(P2-1) P V P2_ ' V^(P2-1)" K V 2/7 

Taking into account equation ( |6.49 ), we arrive at the following result: 



(6.62) 



Proposition 6.2 All constrained Toda lattice eqs.(6.1t)-(6.2l) can be re-expressed as a single set 
of flow equations for one independent matrix Q only: 

d 



Of. 



Q 



>,Q P - 



P2- 1 



8 = l,... ) p 2s 2(p2-l),3(p 2 -l),...,p 1 (p 2 -l) (6.63) 



r(p2-l) 



for r = 2, . . . , pi 



where the "string equation" constraints ( 6.1{\ )-(6.1{> ) become, taking into account (6.56) 







oh dt P2 -i 



'pi— 1 



d 



P2 2 d 



Q , I ^( r + l)t r+1 — + ^ i P2 ^r 1 Q = -l 



r=l 



dt r (P2-1) P2 dh 



(6.64) 



(6.65) 



Similar (dual) equations to d6.63|) -( |6,65 ) can, of course, be obtained for Q as an independent matrix. 

Also, let us write down for the later use the ti-lattice equations of motion for Q = Q^- 1 
resulting from ( |6.63| ) (cf. the analogous Q lattice eqs.( 6.25] )) : 

-b (n) 



d - 1 - ~ 
dh ot P2 _i 



(6.66) 







R,, 



R n (b (n) -b (n - 1)) , S ( n ) = 6 1 (n)-6 1 (n-l) (6.67) 



Similarly, for the higher flows of Q we get the complete analogs of eqs.( 3.38| )-( |5.39 ) : 
J-6 (n) = J- (Q s ) , J-fe!(n) = J- (#») ) 

OT S Oti v y n« OT S ati \ v /n,n+l/ 

L(Q S ) a = l J ...,p 2j 2(p 2 -l),3(p2-l) J ...,pi(p 2 -l 

f-L \ /n+l,n 



(6.68) 



5 7? 
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The last equality in ( |6.66 ) follows from first equation in Q6.68D . 

Comparing the constrained Toda lattice representation ( |6,63| )- (|6.65| ) for the two-matrix model 
with that for the one-matrix model (eqs.(l534D, ([534l) ), we notice their formal resemblance up to 
the following differences: 

• Q in the two-matrix-model case has P2 — 1 > 1 non-zero lower diagonals (this is, of course, 
due to the nonlocal measure in the orthogonal polynomial formalism, cf. (|6.8|)). 



In the two-matrix-model case we have, besides {t r }, an additional subset of fractional flows 

{«.}. 

The fractional flows contribute an additional term to the "string equation" (cf. ( p\34 ) versus 

O). 



6.3 Construction of the Corresponding Continuum KP Hierarchy with the 
String Equation Constraint 

The linear auxiliary problem associated with fl6.63|) -( |6.65 ) reads: 



_id_ 



Mn , T^A 



Q 



P2- 1 

nm 
Pi-1 



A 



dt 



d d 

A = -WrWn 



A( = Qnrn4>m) = (cg-lO)?' 1 + £ PrQ[ 



+) 



P2-1 



A 



(6.69) 
(6.70) 



r=0 



In components, using the parametrization of ( |6.21| ) and ( |6.46| ), we have: 



Mn 
Ot\ 



P2-1 

A+l + a (n)ip n + £ ak(n)ip n _ k 
k=i 



-Rnl/j. 



n-l 



(6.71) 
(6.72) 



Applying again the Bonora-Xiong procedure to ( 6.71| ), (|6.72j ),( |6.7d|) as in derivation of eqs. (|5.20|) — 
(5.21) in the one-matrix-model case, we obtain from ( |6.63| )-( |6,70| ) the following continuum Lax 
problem at a fixed lattice site n (with continuum "space" coordinate x = t\) : 



Mr, 







dt 
_d_ 



T L{n) 



L(n) 



L(n)A , 
L(n) , C s {n) 
<} L(n) 



T^A = -Cs{n)lp n , jr- 

dt s g oX 



A = M(n)i/j n 



, 5=1,... ,P2,2(P2 ~ l),3(p 2 - 1),. • • ,Pl(p2 ~ 1) 
1 



dt P2 -i 



■1 



g 



M(n) , L(n) 



(6.73) 
(6.74) 
(6.75) 



where: 



P2 — 1 



L(n) = a (n) + £ 

k= 

- D^ 1 R n+ \ 



(-l) k a k (n) 



(D x - d x In (h n ■ ■ ■ Rn-k+2)) ■ ■ ■ {p x - d x In R\\ D x 



(6.76) 



47 



s ( _'\\kf)P2- 1 

= a 

fc=l -^n ' ' ' ^Hi-k+l 



L (Ac - da In • • • -R n _ fc+2 ) j ■■■D x 



M(n) 



1 



pi-i 



r=l 



(6.77) 
(6.78) 



Here all coefficients can simply be expressed in terms of the matrix elements R n +i, bo(n), bi(n), . . . 
i 

of Qp?- 1 ( 6.46| ) at a fixed site n through the t\=x lattice equations of motion ( |6,63| ). 

More transparent form for the constrained continuum hierarchy ( 6.74| )~( |6.75| ) is obtained after 
performing a suitable gauge transformation and operator conjugation: 



+ 



b (n) = (L(n))-- 

P2 



L{n) = e/^ o(n) (L{n)X e~^ n) = D^~ x + (p 2 - 1) b x {n)Dl 
+ Rn+tfDx-bo^y 1 
C s {n) = efi o( > n) (Csin))* e~ fi° {n) ° 

M(n) = e/ ?o(n) (M(n))*e-i" ?o(n) 
pi— l 

+ J2 ( r + l)*r+l (L(n)Y {+) 

r=l 



+ ••• + 



PS" 1 + (pi-lfi 2 ( - L ^(+) 1 + 

0(L(n)) 



(6.79) 
(6.80) 

(6.81) 



In derivation of the above equations the expressions ( |6.60| )-( |6.6l| ) were used. The last term 0(L(n)) 
in ( |6.81 ) may be discarded without changing the formalism. 
Correspondingly, ( |6.74| )-( ]6,75| ) acquire the form: 

B_ 

dt s 

d , (> r P2 



' L v v >){+) ' 
8 



n) 



s = 1, . . . , P2 , 2{p 2 - l),3(p2 - 1), • • • ,pi(p2 " 1) (6-82) 



0*1 V ; 9t P2 _i 



a pi- 1 

i)^^ + E ( r + ( L (^))(+) J = 1 ( 6 - 83 ) 



L(n) 



(P2 



The continuum "string equation" constraints ( |6.83| ), upon using the flows from equations ( 6.82] ), 
become (cf. the discrete "string equations" fl6.65| )) : 



OTl d* P2 -i 



'pi— l 



E( r + 1 )*r+i^- + 7^-n^: 







r=l 



dt r (R2-1) "0*1 



(6.84) 



or, in terms of the original evolution parameters using ( 6.6C| ) and eq.( |6.58 ) for s = p 2 — 1 together 
with (ET5C| ) : 

^pi— i q 8 \ 

Y^(r + l)t r+1 —+g-^-\L(n) = -1 (6.85) 

Recalling ( 2.31 ) (or (j2.52|) ), we notice that ( 6.7S| ) is nothing but the Lax operator for the 
cKP P2 _i j i hierarchy which now is subject to the additional "M-constraints" ( 6.84| ). The latter, in 
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turn, are special case of the condition for invariance under the lowest additional symmetry flow 
( 4.451 ) with the order of L being m = p 2 — 1 and the index I = s running over the set of non-zero 



evolution parameters as in eq.( |6.82[) (recall the identification t r = i r ( P2 -i) fl6.64j )). 

We shall denote in the sequel as scKP P2 _i i i the cKP P2 _i 5 i Lax system ( |6.82| ),( |6?f9 ) with the 



additional "string equation" constraint ( |6.83| ). As shown above, it is the exact continuum analog 
of the constrained Toda lattice Lax system fl6.16| )-( |6,20| ) describing the two- matrix model (|6.1|). 

6.4 Partition Function, the Toda Lattice Structure 

Before presenting the explicit formula for the partition function Zn ( |6.1| ) in terms of the r-function 
for the scKPpj-^i hierarchy, let us discuss the role of the Toda lattice site label n in the context of 
the continuum scKP P2 _ii hierarchy. First, we observe that, by construction, the coefficients of the 
Lax operator L(n) ( |6,79| ) satisfy the scKP^.^i system ( |5.82 )-( |6.83 ) for any n = 0, 1, 2 . . . with 



the following "boundary" conditions due to the semi-infiniteness of the underlying Toda matrices 
Q,Q and accounting for their definition ( 6.53| ) : 



Ro = o , M-i) = 7io = - fe " 1)t;2 - 1 



P2t 



in 



731 2 (P2 ~ 2 )t p2 -2 . , 

■z- 7io = J ( 6 - 86 ) 

-3711 P2tp 2 



where the expressions ( 6.6C| )-( |6.61| ) were used. Therefore, Toda lattice shifts n — > n + 1 are equiva- 
lent to mappings of one solution of the scKP P2 _i i hierarchy into another one of the same hierarchy, 
i.e., Toda lattice shifts must correspond to auto-Backlund transformations for the scKP P2 _i.i hier- 
archy. Here we shall show that, indeed, Toda lattice shifts within the discrete constrained integrable 
hierarchy ( |6.63| ) ( |6.65D generate special (constrained) Darboux-Backlund transformations for the 
scKP p2 _i i i system (|6.82;) -(|6^). 

To this end let us rewrite L(n) ( |6.79| ) in the equivalent "eigenfunction" form: 

L(n) = D^- 1 + {p 2 -l)b 1 {n)DP x 2 ~ 3 + --- + ${n+l)D- 1 y{n + l) (6.87) 
$(n + l) = R n+1 exp{Jb (n)} , ^(n + 1) = exp{-y6 (n)} (6.88) 

The Q lattice equations of motion (|6.66j )- (|6,67l ) yield the following recurrence relations for the 
eigenfunctions <&(n), \&(n) ( |6.88 ) and r-function (recall prop.[T^) of scKP P2 _i ) i : 



* d *("+»)_*W_ (6 . 89) 



dhdtp^! *(n) *(n-l) 

f(n+ 1) = ($(n)) _1 -> $(n) = expjy 6 (n)} (6.90) 

J {n) . =exp{ (b (n)X = <5>{n) -» = *(n)*(n - 1) • • • $(0) (6.91) 

T(n — 1) <J J r(— 1) 

Here r(— 1) denotes the r-function of the "initial" Lax operator L(— 1), i.e., L(n) (|6.87| ) for n = — 1, 
which is pure differential due to the "boundary" conditions ( |5.86D : 

L(-l) = Df - 1 - (P2-l)(P2-2)t P2 - 2 £)g2 _3 + . . . (6 Q2) 

V2tp 2 



49 



Lemma 6.1 In terms of the original {t s } flow parameters the "initial" Lax operator L{—1) be- 
comes: 



L(-l) 



P2-1 

£ 

s=0 



s+1 



-t s +l 



d 



t 



P2-1 



dti p 2 t 



(6.93) 



P2 



Proof. The lemma follows easily from ( |6.76j ) and ( |6.79| ) upon using the original "coupling con- 
ditions" ( |6. 19 ) to express the coefficients of L(n) as (using the short-hand notations ( 6.5C| ),( |6"l32|) ) 



a k (n) 



P2-1 



Qn. 



n—k 



Rn ' ' ' Rn—k+1 s=k+l ^ n ' ' ' Rn—k+1 

and subsequently plugging it into (recall x = t%) : 



7 n ct k for n = -1 



(6.94) 



L(-l) 



o 7l0*i 



'P2— 1 
s=0 



(6.95) 



□ 



Comparing ( 6.8S| )-( |6.91| ) with the general formulas of Darboux-Backlund transformations for 
generic cKP r ^ m= i hierarchies ( |3.26| ),( p^36 ),( 3.35| ) allows us to identify constrained Toda lattice site 
shifts with Darboux-Backlund transformations for the scKP p2 _i j i systems as follows: 

L(n) = T(n)L(n - l)T~ l {n) , T{n) = $>(n)D x <5>- l (n) = D x - b Q (n) (6.96) 
$(n + 1) = (T(n)L(n - 1)) $(n) = T(n) • • • T(0) (L(-l) n+1 $(0)) , *(n + 1) = ($(n)) _1 

According to the general formulas ( |3.57| ),( ^59| ), the solutions for the r-function and (adjoint- 

)eigenfunctions of L(n) ( |6,87j ) and, therefore, for the matrix elements of Q = Q^- 1 (cf. ( |6.88j )) are 
expressed entirely in terms of the "initial" eigenfunction ^(O) of the "initial" Lax operator L{— 1) 
^92 



<3?(n) 

T ( n ) 
r(-l) 



w n+1 


$(0),^$(0),. 




w n [$(o), Jr $ (°)>--- 





(6.97) 



n $a) = ^„+i 

d i+ i- 2 <f>(0) 



det 



die 



(6.98) 
(6.99) 



(L(-l))^ 1 $(0) 

S = l,...,P2,2(p2-l) 3 3(p 2 -l),...,Pl(p2-l) 

where we used the constraint ^^(O) = — *^(0) following from the first constraint eq.( 6.84| ). 



Furthermore, the "string equation" constraints (|6.83j ) (or ( |6.84D ) imply the following additional 
conditions for the Lax eigenfunctions 3>(n), which generate the Darboux-Backlund transformations 



'Pi-l 



9 



^(r + l)t r+ i— +ti + £ p2 _i + 



r=l 



L 2 4- I $(n) = 

(P2-1) P2 5t! ' 
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— - = — 



$(n) 



(6.100) 



where the second equality accounts for the identification t T = i r (p 2 -i) Q6-640 . The above equations 
are nothing but a special case (for the current set of evolution parameters ( |6.82[) — ( |6.99[) ) of the 
condition (4.47) for compatibility between additional symmetries and Darboux-Backlund transfor- 
mations of general KP hierarchies. 

To write down more explicit expression for $(0) (|6.99| ) , we shall use the representation (4.50) 
of constrained eigenfunctions in terms of constrained BA functions (satisfying fl4.46| )) which in the 
case under consideration reads (using the short-hand notations ( 6.5C| ), (|6.52|) ) : 



*(0) 

L(-1)^°)(A) 



(6.101) 



' P2 — 1 



AV (0) (A) , i.e. e 710 * 1 £ a slu s D s e" 710 ' 1 V (0) (A) = At/> (0) (A) (6.102) 



d_ 





V (0) (A 



V> (0) (A 



OX 

The solution for ip 



(L(-l))j*^W>(A) , s = 

(L(-l)r^°)(A) = A>(°)(A) 
8 



,P2 



(6.103) 
(6.104) 



8t 



-^ (0) (A) 



P2— 1 

^(r + l)t r+1 — +t 1 +i p2 _ 1 + 



P2 







r=l 



V> (0) (A) 



(6.105) 



°) (A) can be written in the form: 



J^>l A? '*''+ A *P2-l 



7n e 



E r >l Ar *>- +710*1 



,7/i exp |/i (t x + -^-^A ) [ / (0) (/i; {£'}) 
d M e( sA+ ' Y11 * 1 )' i /( )( / x;{t / }) 



(6.106) 



where {t 1 } = fa, . . . ,t P2 ), and we performed change of integration variable p, = 7n/i + 710 to pass 
from the first to the second equality in Q6.106 ). Now, inserting the last expression ( |6,106j ) into 
( 6.102 ) reduces the latter to an equation for f(°\fj,;i') : 



8_ 



/ (0) (M*'}) 



P2 — 1 

(s + i)t s n s + (h - 7 U ti 

s=l 



/ (0) (m;{*'}) 



(6.107) 



Notice, that due to x = ti-independence of 711 = as given in ( |6.60| ), the last term in the 

square brackets does not depend on x = t% . Substituting the solution of eq.( |6.107| ) into ( |6.106 ) we 
obtain: 

(A) = e -«*'})+Wi e E r >i ^ f dfi e g^+E? =1 (6 . 108) 
where the "integration constant" e({i'}) from ( 6.107| ) is determined through the i s -How eqs.( |6.1()^ ) 



at. 



I P2-1 

£ 

(7=1 



<r+l 



Wi) (--*») P2<Ll ^ " " (*i - 711*1 



P2- 1 



(0) 



(6.109) 
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Here the subscript ( ) denotes taking the zero-order term of the corresponding pseudo-differential 
operator which, moreover, has constant x = ti-independent coefficients. Then, for the eigenfunction 
(I6.101D we get: 



*(0) 



M{i'})+~noti 



pi 

I 

-r=l 



1 / d\ l( i\2 exp <^ J2 X l f r + J2 X 2^ + 



(6.110) 



l>2 



where the expression ( |6.8| ) was used. 

In complete analogy to ( |6.11C| ) we can find the relation between h n - the normalization factors 
in the orthogonal polynomial formalism ( |6.8|), a nd the Lax eigenfunctions $(n) in the scKP P2 _i 5 i 
formalism for any n. Namely, consider eq.( 6.90 ) and compare it with the second eq. (|6.23| ) for h n 
which yields (accounting for ( |6.53| )) : 



K = $ 



(n) exp 710 - £ n (i')} 



(6.111) 



where again £ n (t' ) are "integration constants". Employing the lattice equations of motion for h n : 



d 2 



dhdt 



— In h r 



+i 



hr. 



(6.112) 



which follow from ( 6.27|) and ( 6.23[) , fixes the "integration constants" e n {i') = n l n 7n + £ (^') U P to 
an overall n-independent function e(r) which, therefore, must coincide with s(t') from eqs. Q6.10D| )- 
(16.1101) . Thus, we get: 



a>r \ h It (Pa ~ 1 ) t P2-i 
9{n) = h n exp < -t v 



P2t 



e(f) 



P2 



P2 



g{p2 



i) tp2 



(6.113) 



Now, substituting the expressions ( |6.113| ) into eq. Q6.11| ) and taking into account (6.i 
arrive at the final expression for Zn as anticipated in (|6.5D (|6T^) : 



we 



Zjsr 



^-''exp/ivA-^- 1 ^-' 



det 
det 



r(-l) 1 { 
3*+J-2$(0) 



P2t 



s{i>) 



P2 



P2 



g(P2 - 1; 



JV(iV-l) 



(6.114) 



d t+ i- 2 h 



expiNlt ^ 2 ' 1 ^- 1 -e(i') 



P2t 



dti l dti j - 1 



(6.115) 



To get the last equality we used again ( 6.110| ) as well as the chain of simple identities (recall second 
eq.( |6.100| ),( |6.99| ) and the "similarity" form of the initial Lax operator Q6.102| )) : 



c -(tiiio+e(i')) dk 



$(0) 



-(imo+e(f)) (L(-l)) fe $(0) = (e^ l7l0 L(-l)e* l7l0 ) fe / i0 



^9 dh g 



-t\ | h 



0' 



dh 



(6.116) 
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The last equality in ( 6 . 11 5| ) was previously obtained [19] from different approaches. 

Thus, we conclude from ( 6.114 ) ( |6. 1 15| ) that the solution of the discrete two-matrix string 
model ( |6.lD is provided (up to an overall "coupling constant" dependent factor) by the r-function 
of simple DB orbits of the integrable scKP p2 _i^ hierarchy ( |6.79| ),(6.82)-( 6T83"| ) with finite number 
of flows (when p\ is finite). Furthermore, based on arguments parallel to those leading to eq.(|5.4|) 
in the one-matrix model case, we can identify t(N — 1) ( 3.98| ) for generic initial eigenfunctions 
<3?(0), corresponding to more general than ( |6.93| ) initial Lax operators and not subject to "string 
equation" constraint ( 6.100|) , as a partition function of a generalized two-matrix model with matrix 
potentials of non-polynomial type. 



Example: p? = 3 case. Let us illustrate the above general formulas for the simplest nontrivial 
case when the second two-matrix model potential in ( |6.1[ ) is of cubic order, i.e., p2 = 3. The 
corresponding SCKP2,! Lax formulation specializes as follows (recall x = t\) : 



L(n) 



D 2 + 26i(n) + R n+ \ [D-b {n) 

-2 /„- \ -2 



-1 



g / 3t 3 



3\2g J 
2g ( 3t 3 



3*3 
25 



9 



t2 



3*3 
9 



-l 

h 



Eqs.( |6.109| ) specialize to: 







2U 



73- £ (*2»*3) = 



a 



2U 



Of 



3t 3 



01: 



3t 3 



and, therefore, ( 6.110| ) acquires the form: 



ho exp 



r f 2*2 iW \ 



3t 3 6 

r- t 2 2 (t 2 ) 3 

exp<^ t\— — (T 

1 3t 3 3( 3 t 3 ) 2 



3*3 

2ff 



D 2 - 


2t x 

w 3 




h - 


m 


- 2 1 

~9 t2 




f) 


' (•-■ 


*| 




3* 3 



1 (24)' 

5 (*) ! 



Pi 



dX dfi exp < \ r t r + //*s 



,r=l 



8=1 



(6.117) 
(6.118) 

(6.119) 



(6.120) 



(6.121) 



7 Multi-Matrix Models as Constrained KP Hierarchies: Darboux- 
Backlund Solutions 



Similarly to the two-matrix model (|6.1| ), the partition function of the multi- matrix (g-matrix) model 
reads: 

,■ ( <? Pa q-l 

Z N [{t^}, {*^}, {<?}] = JdM!- dM g exp lJ2 E 4? Tr M « a + E ^,a+i Tr M Q M Q+1 



o=l r a =l 



a=l 



(7.1) 

where M a are Hermitian N x N matrices, and the orders of the matrix "potentials" p a may be 
finite or infinite. As in the two-matrix model case, one associates |31] to ( [7.1| ) generalized Toda-like 
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lattice systems subject to specific constraints. Correspondingly, Zn and its derivatives w.r.t. the 
coupling parameters can be expressed in terms of solutions of the underlying Toda-like discrete 
integrable hierarchy where {t^}, • • • , {t^} play the role of "evolution" parameters. 

It turns out that, in order to identify the continuum cKP integrable hierarchy which provides 
the exact solution for ([O]), we need the following subset of the associated linear system and the 



corresponding Lax ("zero-curvature") representation from [31 

d 



d 



dt 







dt\- 



- (Q(i)' 
Q(i)( +)) Q(i 

Q r (+) ,Q(q) 







dt 




(a) 



dti q) 



Q(l)= Q(l),Q(q)i 







at 



Q(q)= Q(q),Q(q): 



(7.2) 
(7.3) 
(7.4) 



Introducing the notations: 
t r = t^> , r = l,...,pi 



f{9) 



1, . . . ,Pg 



Q = Q0) 



Q = Q(q) (7.5) 



the system ([T.2|)- (|7.4[) becomes exactly the same as the one in the two-matrix model case ( |6.14| )- 
( 6.17] ) (with the exception of the second eq.( 6.14| ) which now has a more complicated form but will 
not be needed in the sequel). Furthermore, there is a series of additional constraints ("coupling 
conditions") relating Q = Q(l) and Q = Q(q) which have much more intricate form (involving 
also the "intermediate" Q(2), . . . , Q{q — 1) matrices) than the two-matrix model ones (|6.18| )- (|6.20|) . 
However, their explicit form will not be needed in what follows since we will be able to extract 
the relevant information only from the discrete Lax system ( |7.3|) -([7^) and the relations ( 6.23 ) 
expressing Q = Q(l), Q = Q{q) in terms of orthogonal polynomial factors. 

Here, we shall introduce the same parametrization for the matrix elements of Q = Q(l) and 
Q = Q(q) as in ( |6.21[ )- (|6.22[) (with the obvious changes in the sizes of the Jacobi matrices): 

Qnn = Q(l) nn = Oo , Q(l) n ,„+l = Qn,n+1 = 1 , Q(X)n,n-k = Qn,n~k = a k {n) 



Q(q)nn 
Q(q)n,n+k 



k = l,.. . ,m(l) ; m(l) = (p q - 
„ = Qnm = for m-n>2 , 

n = Qnn = &o( n ) > 0(?)n,n-l 
= Qn,n+k = bk(n)R„^i ■ ■ ■ R 



■1)...(P2-1) 

n — m> m(l) + 1 

= Qn,n—1 — Rn 

-1 
n+k 



(7.6) 



L n+1 

k = 1, . . . , m(q) , m(q) = (p g _! - 1) • • • (p x - 1) 
Q(q) nm = Qnm = for n - m > 2 , m - n > m(q) + 1 (7.7) 

Then, we have the same system of relations ( 6,11 )- (|6.13|) between the Q = Q(l), Q = Q(q) 
matrix elements and the normalization factors in the nonlocally generalized orthogonal polynomial 
formalism [ |60| l (using notations (|7.5| )) : 

h n 5 n 

r=l s=l 
.9-1 Pa 



p(\,H; {t"},{g}) 



{t>>} 



/ dXd^p n (x)e W {Y / \ r tr}p(x^-,{t''},{g}){J2^} p M ( 7 - 8 ) 

J J r=l s=l 

/q— 1 q-1 p a q—2 

- a=2 a =2r a =l a=2 



+ 912^2 + g q -l,qV q -xlx\ 

= (V 2 ),...,^- 1 ) 



(71 
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Furthermore, with the notations ( |7.5| )-( [7?7[) all equations from the two-matrix model case ( |6.25| )~ 
( 6.39| ) are literally satisfied for the (/-matrix model matrices Q = Q(l), Q = Q(q). Also, as in the 
two-matrix case it is natural to introduce the fractional power: 



Q = Q 



m(l) 



Q(l)^ 



(7.10) 



with exactly the same parametrization as in ( 6.46 ) and satisfying the same system of lattice eqs. 
of motion d6.47|) -( |6.48 ). Therefore, we have the following: 

Proposition 7.1 The matrix elements ofQ = Q(l) are completely expressed in terms of the matrix 
elements of Q = Q{q) through the relations: 

m{l) _ . 
Q(-) = E a 4) > Qh =E^H ' « = 0,l,...,m(l) (7.11) 

s=0 cr=0 

where the coefficients ao,7 s ,o are t\ -independent, whereas the coefficients a s ,j sa with a > 1 are 
independent of t\,t\. All -y sa are expressed through a s = 7 m (i) )S exactly in the same way as in the 
two-matrix model case: 



7ss 


= (711 y 


ls,s-2 


= (7n) S ~ 2 


711 


= (a m (i)) 



7s,s-l 

s(s - r 



s(7n) S Sio 



(7io) + s 



f 731 

V.3711 

«m(l)-l 



7io 



(7.12) 
(7.13) 



m(l)-l 

) 



731 

3711 



7io 



m(l) 



a m(l)-2 



m(l)-2 
) 



/ \ m(l 

[«m(l) J 



m(l) (a m (i)) ™ (1 



m(l) - 1 
2m(l) 



) 



/ \ m(l 



(7.14) 



Similarly, we have the exact analog of the dual statement (6.54) with the roles of Q = Q(l) and 
Q — Q(q) interchanged. The difference with the two-matrix case is that now the coefficients a s 
and (3 r do not have any more the simple expressions ( 6.50| ) and ( |6.55| ). 

As an important consequence of (7.11), let us take its diagonal 00-part and use the last eq.( 6.23| ) 
which yields: 



d_ 
dtx 



hn 




+ a \h 



(7.15) 



This equation is the only remnant of the constraints ("coupling conditions") on the multi-matrix 
model Q-matrices which will be used in the sequel. 

Further on we follow the same steps as in the previous section to derive the continuum cKP r m hi- 



erarchy associated with the general (/-matrix model. Namely, as a first step using ( 7.11| ) we re- 



express the Toda-like lattice hierarchy (gj)-(gj) as a single set of flow equations for Q = Q(l) 



m(l) 



3f. 



M>n , T^r^n = ~ (Q(-)) V>m 
Ot s K J nm 

Q,QL , s = 1, . . . ,p q ,2m(l),3m(l), . . . ,p\m(l) 
*rm(l) for r = l,...,p! 



(7.16) 

(7.17) 
(7.18) 
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where, as in the two-matrix case, we have introduced a new subset of evolution parameters |is| 
instead of yt s = ti q ^ defined as: 



d_ 

dt s 



£7. 

(7=1 



_d_ 

Ot n 



,m{q) 



(7.19) 



As a second step, one employs the Bonora-Xiong procedure to get from the discrete Lax system 
(|7.16| )-( 7.17| ) an equivalent continuum Lax system which, upon operator conjugation and similarity 
transformation as in ( |6.79 )— ( p.80|) , acquires the form (as before x = t%) : 

d_ 

-1 



1, 



^L(n)=[(L^(nj) {+) , L(n) 
L(n) = D™V+m(l)b 1 (n)D™V- 2 + -.. + R n+1 (D x -b (n) 



,p ffJ 2m(l),3m(l) > ...,pim(l) (7.20) 



(7.21) 

Exactly as in the two-matrix case, lattice shifts n — > n + 1 in the underlying discrete Toda lattice 
system, described by (|7.2| )-( |7T^) , generate the Darboux-Backlund transformations in the continuum 
cKP, r=m ( 1 ) +l m=1 hierarchy ( 7.20 )-( 7.21 ) and the solutions for the eigenfunctions and r-functions 
at each successive step of the Darboux-Backlund transformation is given explicitly by eqs.( |6.97| )- 
( 6.98| ) where everything is expressed in terms of the eigenfunction of the "initial" Lax operator 
L(— 1). The difference with the two-matrix case is only the explicit form of the latter (cf. ( |6.102j )) 



L(-l) 



O 7ioil 



s=0 



(7.22) 



where the coefficients a a , 710, 711 have more complicated dependence on {t s } than in the two-matrix 
case. 

Exactly as in the two-matrix case, we obtain the relation between the n-th step DB eigenfunction 
$(n) and the orthogonal polynomial normalization factor h n ( |7.8| ) which generalizes ( |6.113| ) : 

$(n) = h n 7i n i exp {t l7 io + e{t')} (7.23) 

Substituting ( |7.23| ) into (6.11) and using the Wronskian formula ( |6.98 ) we get: 



JV-l 



n h - 



det 



n=0 



a 



dti 



/- - \ N(N-l) 

-JV(ti7io+£(*')) 7 ~ 5 



det 







i-i 



i-l 



(7.24) 
(7.25) 



where we absorbed the 711-factors via changing — » by the definition ( |7. 19 ) , i.e., 7n 1 ^" = 



Now, in complete analogy with ( |6.116 ) we find using ( 7.22| ) and ( |7.15| ) : 

'm(l) \ 



(e _fl7l0 L(-l)e* l710 ) J1 ^0 = \J2 



s=0 



9 s 



hr 



-1 



(7.26) 



Substituting ( 7.26 ) into ( [7.25 ) yields the final result for the multi-matrix model partition function: 

d i+ i- 2 h 



Zjsi = det 



(7.27) 



which, as anticipated in ( 6 . 5 ) ( |676|) , is functionally the same as for the two- matrix model ( |6.115| ), 
however, with a more complicated expression for /iq ( |7.8[ ). 
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8 Concluding Remarks 



The intimate interrelation between discrete multi-matrix models and specific Darboux-Backlund or- 
bits of generalized constrained KP integrable hierarchies (cKP rjm hierarchies or, equivalently, SL(r+ 
m,m) KP-KdV hierarchies - see the remark at the end of section |2|) looks now well understood. 
Let us briefly sketch, in conclusion, some directions for further developments connected with the 
present approach. 

First, we will be aimed at extending our construction of the modified additional symmetry 
flows for the cKP r>m hierarchies to cover the whole Wi +oc algebra. The relevant steps were already 
indicated in section Let us recall that eqs.( 4.2 3] ),( |4.27 ) bear direct resemblance to the Wi +oc 



structure, which appeared in |35| where the higher-spin P^-generators have been built up from 
bosonic currents through a generalized Sugawara construction. 

Next, it is desirable to prove a theorem, analogous to the theorem of refs.[42, explicitly 



expressing the action of the modified additional symmetry flows ( 4.29|) (and their generalizations 



for the full Wi +00 algebra) directly in terms of differential operators acting on the r-functions (cf. 

(PD-(PD). 



Recently there was an active interest in multi-matrix models in the context of various random 
matrix problems in condensed matter physics (see refs.|6(| and citations therein). The main object 
of physical interest there is the so called joint distribution function V (Ai, . . . , Aj\r), with the partition 
function being: 



N 

Z N = const / J[ dXj V (Ai, . . . , \ N ) (8.1) 



Recall that linking the discrete (multi-)matrix models with the continuum integrable hierarchies 
(sections |5|47|) amounts to identifying, up to an overall factor, the (multi-)matrix model partition 
functions with T-functions of Darboux-Backlund orbits of cKP rm hierarchies Zjq ~ r^ -1 ) (cf. 
eq,( |6.114l )). In this spirit one associates V (Ai, . . . , Xn) with the integrand in the "eigenvalue" 



(A-integral) representation (3.47). The latter can equivalently be written as: 



J N 

V (Al, . . . , \n) = exp < — -Hone-body (^i) - H two 

— body (-^i i Aj) H man y_ body (Al, . . . , X N ) 

[ J I i>j 

(8.2) 

^onc-body(A) = -ln<rf 0) (A)-$>A' (8.3) 

l>i 

r-l 

^two-body (A,, A,) = -^(Ai-A.-f-ln^AfA^ 1 - 8 ) (8.4) 

s=0 

N 

^many-body(Al,...,A JV ) = - In (t - ]T [Aj 1 ] ) (8.5) 

3=1 

Thus, we see that the generalized multi-matrix models corresponding to the generic DB orbits of 
cKP r i hierarchies, which do not necessarily obey the ordinary "string" equation, yield additional 
new type of contributions to the familiar effective potential in the pertinent joint distribution 
function. These contributions are: the — ln(j)\ (A) term in the one-body potentials ( |8.3| ), the 
second attractive term in the two-body potentials (|3.4D dominating at very long distances over the 



ordinary repulsive first term, as well as an additional genuine many-body potential (8.5) 
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In the same spirit we can interpret the integrand in the A-integral representation ( |3.47| ) of the 
r-function for the most general DB orbit of cKP r]m hierarchies with arbitrary m > 1. It will 
correspond to a joint distribution function of a system of m different types of "particles" . 

It would be interesting to study the physical implications of the emerging new type of joint 
distribution functions (|3.2|)- (|8.5D , especially regarding critical behavior of correlations. On the 
other hand, from mathematical point of view, it is desirable to derive an explicit multi-matrix 
model representation of the r- functions corresponding to generic DB orbits of cKP r m hierarchies 
generalizing the representation ( |5.4| ) for the generic cKPi i T-function. 



A Appendix: Technical Identities. 

We list here for convenience a number of useful technical identities, which have been used extensively 
throughout the text. 

For an arbitrary pseudo-differential operator A we have the following identity: 

xDx-'AxD-'x- 1 ) ^ = xDx^A+xD-'x- 1 - X d x (V^+M) D^x' 1 (A.l) 



where A + is the differential part of A = A + + A_ = Y%=o A i D i + J2-lo ■ 

For a purely differential operator K and arbitrary functions /, g we have the identities: 

(KfD- 1 g)_ = K(f)D- l g , (/ XT 1 g k) _ = f D' 1 K* (g) (A.2) 

Another useful technical identity involves a product of two pseudo-differential operators of the form: 
Xi = fiD- l gi , i = 1,2 : 

X X X 2 = X l {f 2 )D- l g 2 + /!£>-%*(>!) (A.3) 
where: X\(f2) = fid~ 1 (g\f2), etc. . From the above identity it follows the relation |||: 

m k— 1 

L h )_=J2Y, Lk ~ 3 ~ 1 (^)D- 1 (L*) j m (A.4) 

i=l j=0 



for the cKP Lax operator (2.31). 

Let us also list some useful identities involving Darboux-Backlund -like transformation of 
pseudo-differential operators of the Xj-form above: 

T a (QaD^N) T- 1 = ($^V) Z)- 1 ^ 1 (A.5) 

T a (MZ)- 1 ^) 77 1 = MD- 1 (l*(# )) + {T a (Md- 1 (V a <i> a ))}D- 1 <Z>- 1 (A.6) 

T^MD^JvjT; 1 = MD~ X N + {T a (Md- l {N$ a j)}D- 1 <$>- 1 (A.7) 

(z*y^ a ) = -^c 1 (*« w 1 a*.)) (a.8) 

where $ a is one of the eigenfunctions of a cKP Lax operator L ( |2.31| ) and: 

T a = ^D*" 1 , ^>a = K l 

M = T a {M) = $ a d x ($- 1 M) , N^T~ l \N) = -<S>- l d- l {$ a N) 
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Finally, let us recall the following important composition formula for Wronskians |64| : 

W k (f) 



Tk Tk-i 



Ti(f) 



where: 



w 



-D- 



j'-i 



D + In ■ 



Wj- 



j'-i 



W fe _ 1 (/) = W fe [^i,...,Vfe-i ! /] 



(A.9) 

(A.10) 
(A.ll) 
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